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ABSTRACT 


Methods  are  presented  for  calculating  the  laminar  boundary- 
layer  flow  through  separation  to  reattachment  under  the  influence 
of  a  prescribed  pressure  gradient  or,  in  the  case  of  a  supersonic 
main  stream,  under  the  influence  of  "free  interaction"  between  the 
boundary  layer  and  the  main  flow.  In  contrast  to  the  earlier  work 
of  the  authors  using  the  Tani  method  and  quartic  profiles,  the 
present  method  is  baseo  on  the  Dorodnitsyn  method  of  integral 
relations  and  uses  a  rational  velocity  profile  which  accounts 
properly  for  the  separation  singularity.  As  a  result,  the  possi¬ 
bility  of  higher  approximations  is  inherent  in  the  method.  The 
calculated  solution  for  free  interaction  goes  smoothly  through 
the  separation  point  and  is  in  good  agreement  with  certain  fea¬ 
tures  of  the  Navier-Stokes  solution  in  the  neighborhood  of  separa¬ 
tion.  Good  agreement  is  exhibited  between  experimental  and  cal¬ 
culated  pressure  distributions  up  to  reattachment  for  the  several 
cases  for  which  the  comparisons  were  made.  A  computer  program 
based  on  the  work  has  been  prepared  for  two-dimensional  flow.  It 
is  planned  to  continue  the  work  to  cover  nonadiabatic  boundary 
layers  and  axj symmetric  bodies. 
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CALCULATION  OF  LAMINAR  SEPARATION  WITH 
FREE  INTERACTION  BY  THE  METHOD  OF 
INTEGRAL  RELATIONS 

PART  I  -  TWO-DIMENSIONAL  SUPERSONIC  ADIABATIC  FLOW 

1 .  INTRODUCTION 

This  report  deals  with  automatic  computation  studies  which 
have  been  directed  toward  the  calculation  of  laminar  separated 
boundary-layer  flows  at  subsonic  and  supersonic  speeds.  The 
objective  of  the  present  study  is  the  development  of  a  computer 
program  for  determining  the  viscous  field  including  the  separa¬ 
tion  and  reattachment  zones  over  flat  plates  with  compression 
surfaces.  In  Reference  1,  Abbott,  Holt,  and  Nielsen  developed  a 
method  for  calculating  such  a  flow  field  using  the  method  of  Tani 
and  utilizing  quartic  velocity  profiles.  While  the  method  pre¬ 
dicted  the  experimental  pressure  distribution  accurately  up  to 
separation  and  slightly  beyond,  it  had  definite  limitations  with 
regard  to  the  region  downstream  of  separation  and  with  regard  to 
cooled  boundary  layers.  In  Reference  2,  Lees  and  Reeves  used  the 
method  of  Tani,  together  with  similarity  velocity  profiles.  In 
this  report  the  authors  abandon  entirely  the  method  of  Tani  and 
adopt  the  Dorodnitsyn  method  of  integral  relations  of  Reference  4. 

The  general  approach  of  the  foregoing  methods  is  to  make  some 
plausible  assumption  regarding  velocity  profiles  which  is  sufficient 
to  reduce  the  integro-dif ferential  equations  of  the  boundary  layer 
to  ordinary  differential  equations  which  can  be  readily  solved 
numerically.  For  instance,  in  Reference  1,  a  one-parameter  family 
of  quartic  profiles  was  used  while  in  Reference  2,  a  one-parameter 
family  of  Stewartson  reverse  flow  profiles  was  used.  It  seems 
that  the  critical  factor  determining  the  degree  of  success  enjoyed 
by  any  one-parameter  method  for  any  particular  case  is  the  degree 
to  which  the  assumed  velocity  profiles  will  fit  the  nonsimilar 
family  of  profiles  which  actually  occurs. 

However,  a  shortcoming  of  a  one-parameter  family  of  profiles 
in  the  existing  methods  is  that  they  do  not  permit  higher 
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approximations  to  the  solution  and  thus  do  not  permit  assessment 
of  how  closely  the  assumed  velocity  profiles  approximate  the  actual 
nonsimilar  family  of  profiles.  Another  shortcoming  is  that  one 
cannot  expect  a  one-parameter  family  of  velocity  profiles  to  repre¬ 
sent  accurately  all  possible  velocity  profiles  that  can  be  developed 
in  separated  and  attached  flows. 

In  accordance  with  these  critical  remarks  concerning  one- 
parameter  families  of  profiles,  the  method  of  the  previous  investi¬ 
gation  (Ref.  1)  has  been  completely  abandoned  and  the  Dorodnitsyn 
method  of  integral  relations  (Ref.  4)  has  been  adopted.  This 
method  contains  within  its  mathematical  framework  the  possibility 
of  carrying  out  higher  approximations,  and  with  the  proper  treat¬ 
ment  of  the  singularities  appears  to  contain  the  possibility  of 
approximating  exact  solutions  as  closely  as  possible.  In  the 
present  work  a  fourth  approximation  is  used  so  that  the  velocity 
profile  used  in  a  four-parameter  profile. 

Among  the  questions  to  be  discussed  is  the  interesting  one 
of  whether  the  boundary-layer  equations  give  an  adequate  represen¬ 
tation  of  the  flow  in  the  neighborhood  of  the  separation  point, 
or  whether  the  Navier-Stokes  equations  are  required.  Also,  some 
light  is  shed  on  the  question  whether  the  concept  of  free  inter¬ 
action  can  be  expected  to  apply  with  any  degree  of  precision  in 
the  neighborhood  of  the  separation  and  reattachment  points.  "Free 
interaction"  implies  that  the  body  contour  is  increased  by  the 
boundary-layer  displacement  thickness  when  using  inviscid  super¬ 
sonic  flow  theory  to  calculate  the  pressure  distribution  acting 
on  the  boundary  layer. 
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2 .  GENERAL  CONSIDERATIONS 
2 . 1  Description  of  Problem 

The  problem  under  consideration  is  that  of  calculating 
accurately  by  automatic  means  laminar  separated  flows  under  the 
action  of  prescribed  pressure  gradients  or  under  the  influence  of 
free  interaction  between  the  boundary  layer  and  the  supersonic 
outer  flow.  By  "free  interaction"  we  consider  that  the  pressure 
distribution  of  the  outer  flow  is  the  result  of  a  mutual  inter¬ 
action  between  the  boundary  layer  and  the  outer  flow.  The  bound¬ 
ary  layer,  on  the  one  hand,  through  its  displacement  thickness 
turns  the  outer  flow;  and  the  outer  flow,  on  the  other  hand, 
through  its  pressure  rise  tends  to  thicken  and  separate  the 
boundary  layer. 

Usually  the  boundary-layer  equations  are  solved  in  a  step- 
by-step  integration  as  a  straightforward  parabolic  boundary-value 
problem,  despite  the  fact  that  pressure  waves  in  the  subsonic  part 
of  the  boundary  layer  always  make  the  boundary-value  problem  ellip¬ 
tic.  In  a  sense  it  is  a  corollary  of  the  free-interaction  concept 
that  the  separated  flow  is  independent  of  the  direct  influence  of 
the  downstream  configuration.  By  this  we  mean  that  the  downstream 
configuration  cannot  influence  steady  separated  flow  by  means  of 
upstream  pressure  waves,  although  it  is  thought  that  the  downstream 
configuration  actually  fixes  the  position  of  separation  in  accord¬ 
ance  with  some  reattachment  criterion.  The  correct  position  of 
the  separation  point  must  be  found  by  iteration  in  accordance 
with  a  reattachment  criterion,  and  in  this  sense  the  boundary- 
value  problem  is  elliptic.  However,  once  the  separation  point 
has  been  specified,  then  solution  of  the  boundary-layer  equations 
including  free-interaction  is  independent  of  the  direct  influence 
of  the  downstream  configuration  until  the  calculational  marching 
procedure  reaches  the  configuration. 

The  particular  configuration  which  we  will  consider  is  shown 
in  Figure  1.  Herein  a  ramp  induces  separation  on  a  flat  plate 
and  reattachment  occurs  on  the  ramp.  This  paper  will  present 
results  for  the  pre-separation,  post-separation,  and  reattachment 
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regions.  The  flow  beneath  the  u  =  0  line  will  be  called  the 
inner  flow,  and  that  above  the  u  =  0  line,  the  outer  flow. 

With  regard  to  the  slopes  of  the  compression  surface,  the 
assumption  has  been  made  that  these  are  small  so  that  the  boundary- 
layer  equations  in  a  coordinate  system  parallel  and  perpendicular 
to  the  flat  plate  can  be  retained  all  the  way  to  reattachment  (see 
Fig.  1).  The  additional  terms  introduced  into  the  boundary-layer 
equations  in  the  reattach^nent  region  by  not  taking  the  vertical 
axis  normal  to  the  ramp  are  proportional  to  the  slope  of  the  ramp 
and  are  therefore  negligible  for  small  enough  slopes. 

2.2  Assumptions 

A  number  of  assumptions  are  made  for  the  purpose  of  this 
paper. 

(1)  The  governing  equations  are  those  for  a  compressible 
laminar  boundary  layer. 

(2)  The  Prandtl  number  is  unity. 

(3)  The  air  behaves  as  an  ideal  gas. 

(4)  The  wall  is  adiabatic  and  at  stagnation  temperature. 

(5)  Within  the  range  of  interest  the  viscosity  varies 
linearly  with  temperature. 

(6)  The  pressure  at  the  outer  edge  of  the  boundary  layer 
is  governed  by  the  Prandtl -Meyer  relationship. 

(7)  There  is  no  upstream  influence  due  to  pressure  waves 
in  the  boundary  layer. 

(8)  A  compression  ramp  with  small  slopes  but  otherwise  of 
arbitrary  shape  induces  separation  of  the  boundary  layer. 

2.3  Partial  Differential  Equations  and  Boundary  Conditions 
2.3.1  Physical  plane 

The  laminar  boundary-layer  equations  being  used  are  the 
usual  ones. 


pu 


dx 


+  pv 


du  _  _  .&E  . 

dy  ox 


(1) 
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hmi  +  ii£vi .  o 

ox  oy 


The  domain  within  which  a  solution  to  the  boundary-layer 
equations  is  being  sought  is  shown  in  Figure  2. 

Boundary  A:  w(x)  =  0;  xq  <*.  x  xc  (3) 

w(x)  =  given  function;  xc  <  x  x^. 

Boundary  B:  x=xq  0  <  y  ® 

Boundary  C:  y  =  °°  £  *  £  xr 

The  velocity  boundary  conditions  are: 


Boundary  A: 

u  -  0  v 

=  0 

(4) 

Boundary  B: 

u  =  ui(y) 

u  =  u  at 

o 

y  =  oo 

(5) 

0  at 
dy 

y  =  oo 

Boundary  C:  (a)  Prescribed  pressure  distribution: 

u  -  (x)  (6) 

u.  =  u  at  x  =  x 
i  o  o 


(b)  Free  interaction: 

u  =  u  at  x  =  x  (7) 

1  o  o 

The  pressure  boundary  conditions  include  the  usual  boundary- 
layer  assumption  that  p  depends  only  on  x. 

For  prescribed  pressure  distribution: 

p  ■  p  (x) ,  a  given  function  (8) 
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For  free  interaction: 


P  ■ 


at 


x  ” 


x 

o 


(9) 


The  total  temperature  is  constant  throughout  the  flow 


T  (l  ♦  M") 


(10) 


2.3.2  Stewartson  plane 


In  introducing  the  Stewartson  transformation,  Reference  3, 
we  first  make  a  nonlinear  stretching  in  x  and  y  to  X  and  Y 
as  follows: 


X 


o 


Piai 

poao 


dx 


Y 


Piai 

PoaO 


(ID 


Under  this  transformation,  it  can  be  shown  that  Equations  (1) 
and  (2)  transform  to 


where 


U 


+  v  u 

dX  SY  i 


+ 


2U  £v 

dx  dY 


0 


(12) 

(13) 


(14a) 


V 


Piax 

P  a 

ro  o 


P  a 
p  a 

i 


(14b) 


It  is  noted  that  the  U  profile  is  simply  a  uniform  stretch¬ 
ing  or  shortening  of  the  u  profile  for  a  fixed  value  of  x.  The 
V  profile  is  a  nonlinear  stretching  or  shortening  of  the  v  pro¬ 
file  for  fixed  x. 
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Notice  that  the  entire  Boundary  A  goes  into  Y  =  0  and 
Boundary  C  goes  into  Y  =  00 .  Boundary  B  goes  into  X  =  Xq.  We 
thus  have  the  velocity  boundary  conditions: 


Boundary  A:  Y  =  0 

Boundary  B :  X  =  Xq 

Boundary  C:  Y  =  °° 


U  =  0  V  =  0 

U  =  U^Y) 

U  '  °1  '  “l  if  ^  -  0 


(15) 

(16) 

(17) 


?.3.3  Dorodnitsyn  plane 

Following  Dorodnitsyn  the  physical  coordinates  £  and  t) 
are  introduced  as  follows: 


o 


^ =  ri 

o 


(18) 


New  velocity  variables  are  introduced  as  follows 


U  ~  _  V_ 

u  ■  u  ’  v  '  U, 


(19) 


The  partial  differential  equations  then  become 


at +  ^  V  +  UTi  u7  >"  0 


u. 


-  &L  X  f  “  4.  “  1 

u  di  +lv  +  UT1 0“ 


U  n2- 

&  =  (1  -  u2)  +  *L_u 

aTi  ui  an2 


>  (20) 


Introducing  the  new  variable 


-  -  -  Ui 

W  =  V  +  UT)  — 


(21) 


n 
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we  have 


0 


-  du  ,  —  du 

u  af  +  w  arj 


(1  -  u2) 


(22a) 

(22b) 


The  dot  denotes  differentiation  with  respect  to  £ . 

Under  the  Dorodnitsyn  transformation,  the  velocity  boundary 
conditions  are  transformed  as  follows: 


Boundary  A:  t\  =  0 ;  u  =  v  =  w  =  0 

Boundary  B:  £  =  £  ;  u  =  u^(r]) 

Boundary  C:  rj  =  °° ;  u  =  1;  =  0 


(23) 
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3.  REDUCTION  OF  INTEGRO-DIFFERENTIAL  EQUATIONS  TO  ORDINARY 
DIFFERENTIAL  EQUATIONS 

3.1  Derivation  of  Integro-Dif f erential  Equations 

The  scheme  used  to  obtain  solutions  to  the  foregoing  boundary- 
value  problem  is  to  derive  from  the  system  of  partial  differential 
equations  given  by  Equation  (22)  a  set  of  ordinary  differential 
equations,  the  solution  of  which  will  be  a  close  approximation  to 
the  exact  solution.  The  possibility  of  increasing  the  accuracy 
of  the  approximation  is  also  a  feature  specially  being  sought. 

The  method  of  integral  relations  of  Dorodnitsyn,  Reference  4, 
seems  to  be  one  containing  these  features.  For  this  purpose  we 
derive  the  integro-dif ferential  equations  for  the  inner  and  outer 
layer  from  which  the  ordinary  differential  equations  are  derived. 

Let  us  first  introduce  a  boundary  s(x)  in  the  physical 
plane,  Figure  1,  which  connects  the  separation  point  and  reattach¬ 
ment  point.  Let  u  be  zero  along  this  line  so  that  the  direct 
flow  and  the  reverse  flow  are  separated  into  two  regions  called 
the  outer  and  inner  regions,  respectively.  Consider  first  the 
outer  region  in  the  Dorodnitsyn  plane  in  which  the  boundary  s (x) 

transforms  into  r)  (^).  Introduce  a  family  of  smoothing  functions 
_  s 

f  (u)  given  as  follows: 
n 

fn(u)  =  (1  -  u)n  (24) 


Multiply  the  continuity  equation  (Eq.  (22a))  by  f(u)  and  the 
momentum  equation  (Eq.  (22b))  by  f'(u),  add  and  integrate  between 

n  and  00  to  obtain 

's 


00  00 

”>8  ''s 


00  u  00 

J  f*  jp  (1  -  u2)dr|  +  J  f 

^s  1  ^s 


a2u 

dq2 


dq 


(25) 


Let  us  now  consider  the  velocity  boundary  conditions  on  the  line 

s (x)  as  it  transforms  into  Y  (X)  in  the  Stewartson  plane  and 

s 

t)  (£)  in  the  Dorodnitsyn  plane 
s 
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s  (x)  : 


u  -  0 


VX)* 


U 


0 
p  a 

*  rv 


V 


o  o  P(x.s) 


P  a 

l 

s  (x) 


/ 


P  a  n 

k  JL 


w  (x) 


p  u 
1  o  o 


p. 


dy 


% (£)  :  u  =  0 


W  =  v 
s  s 

Ui 

r‘s  =  uT  V  Ys 
o  V  o 

We  can  now  show  easily  that 


(26) 


I  (27) 


(28) 


CX)  00 

/  st  (£:r)dT1  =  dV  /  £“  dT) 

1s  I* 


/  ^  (fw)d1 


-f  v  since  f (» ) 
s  s 


0 


>  (29) 


f 

J 

% 


f ' 


>2 


U 


dif 


dii 


f. 

s  dii 


j 


Equation  (25)  thus  becomes 
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J  fudrj-  fsvs 


{• 


(30) 


This  is  the  integro-dif f erentia]  equation  of  the  outer  flow. 

With  regard  to  inner  flow,  let  us  integrate  between  q  =  0 
and  Tj  =  i)  to  obtain 


s  U. 


J  ^  (fu)dt,  +  J  (fw)dTl  -  J  V  ^ 


—  (1  -  u  )dq 


'I 

/ 


+  /  f '  ^-7  dq  (31) 

di] 


It  is  then  easy  to  see  that 


^  f"  fu  dr,  +  £svs  -  fS  f  ^  (1  -  u2)dn 


s  s 


+  f  t* 

s  dq 


_  f. 

w  dq 


(32) 


This  is  the  integro-dif ferential  equation  of  the  inner  flow. 

If  we  let  du/dq  be  a  function  only  of  u,  Equations  (30) 
and  (32)  become 


_d_  f  fu  du 

d4  J  £ll 

us  a i) 


s  dq 


1 
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fu  du 
hi) 


+  f  V 
s  s 


(1  -_u‘  )  du  +  f  , 
du  s 


di) 


u 

s 


(34) 


Note  that  u  is  zero  but  is  retained  in  its  notational  form  to 
s  _ 

indicate  the  integration  of  Equation  (33)  is  from  the  u  =  0  line 
in  the  flow  to  u  =  1,  while  the  integration  of  Equation  (34)  is 
between  two  lines  on  which  u  =  0. 


3.2  Velocity  Profiles 

The  choice  of  velocity  profiles  is  a  critical  item  in  the 
success  of  the  present  method.  In  Reference  4,  Dorodnitsyn  uses 
a  representation  for  the  velocity  gradient  du/dT)  in  terms  of 
u  in  such  a  fashion  that  if  r]  is  zero,  the  integrals  in  Equa- 
tion  (30)  are  all  converted  to  integrals  of  u  with  limits  of 
0  and  1.  Consider  only  the  outer  flow  and  note  the  following: 

At  t)  =  00 :  u  =  1  du/dr]  =  0 

(35) 

At  separation,  tj  =  0:  u  =  0  du/dr)  -  0 

If  we  assume  that  the  zero  at  u  =  1  is  a  simple  zero,  neither 
integral  with  du/dr,  in  the  denominator  will  be  divergent  because 
the  particular  choice  of  smoothing  function  has  a  compensating 
zero  in  the  numerator.  However,  if  the  zero  at  separation  is  a 
simple  zero,  then  the  integrals  will  be  divergent.  If  we  choose 
the  zero  in  du/d^  at  separation  as  a  square  root,  the  integrals 
will  be  convergent.  Also,  this  form  is  compatible  with  the  fact 
that  du/dr)  is  a  double- .  alued  function  of  u  just  downstream 
of  separation.  Accordingly,  we  choose  the  following  form  for  the 
velocity  profile  for  the  outer  flow: 

(1  -  u )-yfu  +  c 

- - - - -  (36) 

c  +  c  u  +  c  ,u 
o  1 


dn 
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In  this  formulation  c  ,  c, ,  c  ,  c  ,  • • •  are  taken  as  functions 
of  £  only.  Equation  (36)  is  introduced  in  Equation  (33)  and  as 
many  of  the  f^tu)  smoothing  functions  are  used  as  there  are  c’s. 
Since  the  known  singularities  are  explicitly  exhibited  and  the 
"smooth  part"  is  represented  by  a  Taylor  series,  the  accuracy  of 
the  approximation  should  increase  as  the  number  of  terms  in  the 
series  is  increased. 

With  regard  to  the  inner  flow,  there  is  a  zero  of  du/dr)  at 
the  point  of  maximum  velocity  in  the  reverse  flow.  This  zero 
causes  divergent  integrals  in  Equation  (34) .  Accordingly,  Equa¬ 
tion  (32)  will  be  used  for  the  inner  flow  since  it  will  not  be 
divergent  if  u  is  expanded  in  a  power  series  in  rj  in  the  usual 
way.  Let  us  consider  the  following  boundary  conditions  for  the 
inner  flow: 


n  -  0 

n  =  t,s 


u  -  0 

u  =  0 


hn 

dT] 

dii 


-  a 


w 


a 


(37) 


The  velocity  can  then  be  written  as  a  cubic  equation 


u  =  a  n 
w  1 


(2a  +  a  ) 

- * - 8-  n2  + 

^s 


(a  +  a  ) 
w  s  S 

- r-  v 


(38) 


3.3  Integral  Relations  for  Outer  Flow 

In  the  following  work  we  will  use  a  fourth  approximation  for 
the  outer  flow  so  that  four  ordinary  differential  equations  for 
c  ,  ci ,  c2,  and  c3  will  be  required.  Through  the  use  of  Equa¬ 
tions  (24),  (33),  and  (36),  we  obtain 
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ik. 


d  /  (1  -  u) n  'L  u 


+  c 


(c  +cu+cu2)du-v 
\  O  1  2  J  s 


Ui  j  n  ( 1  -  u)  n_1  (1  -f  u) 

u  . *  . . . . 


V* 


+  c. 


c  +  c  u  +  c„u2  )  du 
O  1  2  / 


-  n,/=7 

ii  —  —  * 


o 

1 


n(n  -  1)  (1  -  u)n-i 


+  c 


-  du 


(39) 


C  +  c  u  +  C0U 
O  1  2 


It  is  to  be  noted  that  all  the  integrals  are  tractable.  If  a 
fifth  approximation  were  used,  the  last  integral  would  contain  a 
cubic  in  the  denominator,  and  the  definite  integrals  would  involve 
elliptic  integrals.  If  higher  approximations  than  the  fourth  are 
desired,  it  is  probably  a  practical  necessity  to  evaluate  the 
integral  by  numerical  means  rather  than  analytical  means. 

To  obtain  vg  we  make  use  of  the  inner  flow  velocity  pro¬ 
file.  From  continuity  we  can  write 


tr 


I 

r 

i 


(40) 


Substituting  u  from  Equation  (38)  yields  for  v^ 

\  ~  12  d\  [V(C,S  -  aw>] 

wherein  from  Equations  (36)  and  (37) 


(41) 


(42) 


since  du/d'')  is  to  be  continuous  across  the  u  =  0  line. 

Let  us  define  two  families  of  definite  integrals  encountered 
in  Equation  (39)  by  the  following  notation: 
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9n(c3>  = 


— n 

u 


V 


du 


(43) 


u  +  c. 


P  (c  ,  c ,  ,  c  ,c  )  = 
n  o’  i ’  2  *  ^ 


u1"1 -yj~u  +  cj  du 


c  +  c  ,  u  +  c  u 
O  1  2 


(44) 


Formulas  for  evaluating  g  (c  )  are  derived  in  Appendix  I.  The 

n  >3 

integrals  are  evaluated  numerically  because  of  the  numerous 


subcases  in  the  analytical  integration.  A  set  of  four  equations 
for  cq>  cx>  c2,  and  c3  is  obtained  from  Equation  (39)  for  the 
four  smoothing  functions  with  n  =  1,  2,  3,  and  4.  This  set  can 
be  somewhat  simplified  by  addition  or  subtraction  of  the  four 
equations  to  yield  the  following  set: 

9^0  +  g2=i  +  93c2  +  (V,  +  +  c2g3)i3  +  [cogo  +  (cQ  +  c.Jg, 

n  u  -\Tc~ 

+  (c,  +  c2)g2  -  c2g3J  jj-  -  v  -  — 

J  1  o 

(45) 

g„c  +gc  +  g  +  (c  g„  +  c  g  +  c„q  )c  -[eg 

^2  0  ^3  1  ^4  2  o^2  1*3  2  J  4  3  [_  0*0 

U 

+  (C1  '  C0)91  +  (C2  ‘  C1  '  2tO)g2  "  (2C1  +  C=,93  -  2C29-J  U7 

*  -  c  +  2(p0  -  V 


(46) 


g3co  +  g4cx  +  g5c2  +  (cog3  +  cxg4  +  c2g5)c3  +  ["2cogi 

u 

-  (2ox  -  co)g2  +  (3c0  +  Cx  -  2c2)g3  +  (3cx  +  c2)g„  +  3c2g5j  — 


-2P  +  8P  -  6P„ 

012 
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c  +  g  c  +  q  c  +  (c  g  +cg  +  c  g  )c  -  |+3c  q 

40  ^51  2  0^4  iy5  2^e'  3  |_  rr*  2 

(-co  4-  3C;L)g3  +  (3c2  -  cx  -  4co)g4  -  (c2  +  4Cl)gb  -  4cpgeJ  — 


=  -6?^^  +  18P2  -  12P3 

(48) 


In  Equation  (45)  the  value  of  v  evaluated  using  Equation 

s 


(41)  is 


v 

s 


t)  %% 


(49) 


It  is  noted  chat  Equations  (45)  to  (49)  are  a  simultaneous  set  of 
linear  ordinary  differential  equations  for  c^,  cx ,  c  ,  c3,  and 
Ui .  For  a  prescribed  pressure  distribution  U1  is  known,  and 
Equations  (45)  to  (48)  are  adequate  for  evaluating  the  boundary- 
layer  flow  up  to  separation.  When  free  interaction  is  assumed, 
an  additional  equation  is  required. 

3.4  Integral  Relation  for  the  Inner  Flow 

Equation  (32)  for  the  inner  layer  involves  a  finite  range  of 
integration  and  contains  no  singularities.  No  particular  advantage 
occurs,  therefore,  in  transforming  to  u  from  t).  In  fact 
divergent  integrals  are  so  introduced.  Furthermore,  there  is  no 
necessity  to  use  the  same  smoothing  functions  used  for  the  outer 
flow.  However,  it  is  convenient  to  use  1  -  u  as  the  smoothing 
function  in  Equation  (32),  which  now  becomes 


u) u  dn  +  v 

s 


(a  -  a  )  (50) 
s  w 


Through  Equation  (40)  we  further  simplify  to 
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u 


dTj  = 


(1  -  u2)dr]  -  (ag  -  aw) 


(51) 


Carrying  out  the  integrations  yields  the  inner  flow  equation 


(52) 


3.5  Free-Interaction  Equation 

The  free-interaction  equation  is  one  relating  the  axial  pres¬ 
sure  gradient  to  the  rate  of  increase  of  the  boundary-layer  dis¬ 
placement  thickness.  The  increase  in  d5*/dx  turns  the  external 
flow  causing  a  pressure  rise.  The  pressure  rise  tends  to  cause 
further  increases  in  5*.  We  will  assume  that  the  pressure  dis¬ 
tribution  is  related  to  the  turning  angle  <p  of  the  external 
flow  by  the  Prandtl-Meyer  relationship  for  inviscid  compressible 
two-dimensional  flow. 


Pi 


•yM  2  dtp 


cos 


-  COS  0 


(53) 


The  turning  angle  0  is  determined  by  the  slope  dw/dx  of  the 
boundary  and  the  slopes  d6Q*/dx  and  d5^*/clx  due  to  the  dis¬ 
placement  thicknesr  of  the  outer  and  inner  flows.  The  relation¬ 
ship  is 
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db  * 
o 

dx 


db.  * 


+ 


1 

dx 


tan  0 


(54) 


It  is  Equations  (53)  and  (54)  together  which  govern  free  interaction 
between  the  outer  flow  and  the  boundary  layer.  The  values  of  the 
boundary-layer  terms  in  Equation  (54)  are  derived  in  Appendix  II 
for  the  assumed  profiles. 
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4.  PRE-SEPARATION  ANALYSIS,  PRESCRIBED  PRESSURE  DISTRIBUTION 

In  order  to  see  how  accurately  the  present  method  will  predict 
nonsimilar  velocity  profiles  in  the  face  of  an  adverse  pressure 
gradient,  it  would  be  desirable  to  compare  it  with  some  solution 
of  known  accuracy.  A  suitable  example  for  comparative  purposes 
is  given  in  the  work  of  Hartree  and  of  Leigh,  References  5  and  6. 
These  particular  numerical  solutions  apply  to  an  incompressible 
laminar  boundary  layer  in  which  the  main  stream  produces  a  pressure 
gradient  with  a  free-stream  velocity  given  by 

Ui(x)  x 

U  (X)  1  "  8l 

o 

We  will  now  apply  Equations  (45)  to  (48)  to  the  solution  of  this 
problem,  and  then  make  a  comparison  with  the  previous  numerical 
solutions.  The  proper  initial  conditions  for  this  solution  are 
the  values  of  cQ,  cx ,  c^,  and  c3  which  best  represent  the 
Blasius  profile. 

4.1  Initial  Conditions  for  Velocity  Profile 

It  is  necessary  to  specify  an  initial  velocity  profile  in  the 
6, ,  r)  plane.  The  initial  conditions  are  specified  at  some  point 
in  front  of  the  separation  point,  and  for  i  <  £o,  it  is 
assumed  that  the  static  pressure  is  constant.  In  a  real  flow, 
interaction  starts  at  the  leading  edge,  however,  we  will  not 
start  the  solution  at  ^  =  0  in  order  to  avoid  the  leading-edge 
problem  associated  with  infinite  shear.  The  initial  velocity 
profile  in  the  £,,rj  plane  under  the  assumed  condition  will  be 
the  Blasius  profile.  In  cases  where  induced  pressure  gradients 
do  occur  at  the  leading  edge  of  a  plate,  the  effects  of  the  pres¬ 
sure  gradients  on  the  velocity  profile  will  diminish  as  the  down¬ 
stream  distance  increases,  anc  the  profile  will  approach  the 
Blasius  profile. 

Let  us  consider  the  analog  of  the  Blasius  solution  in  the 
Dorodnitsyn  plane.  The  equations  in  the  X,Y  plane  are  those 
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for  incompressible  laminar  flow.  For  such  flow  with  no  pressure 
gradient,  we  have 


U  =  U  constant 
o  1 


p  =  p  constant 
ro  r  1 


With  these  relationships  the  Blasius  profile  can  now  be  developed 
using  Equations  (45)  to  (48)  provided  a  suitable  set  of  initial 
conditions  are  found  for  these  equations.  However,  the  suitable 
initial  profile  is  the  Blasius  profile  itself.  We  overcome  this 
dilemma  by  noting  that  a  laminar  boundary  layer  on  a  flat  plate 
with  no  pressure  gradient  will  develop  the  Bias '.us  profile  asymp¬ 
totically  even  if  the  initial  conditions  vary  substantially  from 
the  Blasius  profile.  Accordingly,  a  rough  approximation  to  the 
values  of  cq,  c  ,  c2 ,  and  c^  for  the  Blasius  profile  was  used 
for  the  initial  conditions,  and  the  true  Blasius  profile  was 
sought  asymptotically. 

The  question  arises  how  the  asymptotic  profile  will  be 
recognized  in  the  computed  results.  We  note  that  the  similarity 
of  the  Blasius  profile  states  that,  for  constant  Reynolds  number 
in  the  free  stream,  the  value  of  Y  for  a  given  value  of  u 
increases  as  X  .  ^or  the  slope  du/dr],  we  can  thus  write 
asymptotically 
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The  only  way  that  the  desired  result  can  be  approached  is  if 


c 


(57) 


as  4  — *■  00  • 

In  Table  I  the  initial  values  of  c  ,  c.  ,  c  ,  and  c  are 
shown  and  the  above  ratios  are  tabulated  as  a  function  of  I-  to 
see  if  a  similar  solution  is  developing.  It  is  seen  that  a  simi¬ 
lar  solution  does  develop  vith  the  following  approximate  values 

k  =  3.157,  k,  =  -1.923,  kP  =  -0.3133,  c  =  1.1000 

It  now  remains  to  compare  the  profile  obtained  this  way  with 
the  Blasius  profile.  For  this  purpose  we  solve  Equation  (36), 


k  +  k  u  +  k„u2 
_2 _ 1  .  2..— 

(1  -  u)i/u  +  c3 


du 


(58) 


In  this  way  the  value  of  q  for  a  given  u  can  be  readily  cal- 

i/2 

culated.  The  velocity  gradients  vary  as  £  so  that  we  have 


(1  -  u )  u  +  c3 

k  +  k  u  +  k  u2 
o  i  2 


(59) 


The  values  of  the  velocity  profile  and  its  derivative  as  cal¬ 
culated  by  the  preceding  equations  are  tabulated  together  with 
those  of  Blasius  in  Table  II.  The  Blasius  values  are  from  Refer¬ 
ence  7,  page  107.  From  the  table  it  can  be  seen  that  for  a  given 
value  of  Y,  the  values  of  u  do  not  differ  by  more  than  0.002. 
This  accuracy  is  considered  adequate  for  specifying  the  initial 
velocity  profile. 
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The  skin  friction  calculated  from  the  present  results  is  a 
close  approximation  to  that  obtained  by  the  Blasius  formula.  We 
thus  have 


T 


0.332  p  U  2 
ro  o 


(60) 


4.2  Calculative  Example 

Utilizing  the  initial  values  shown  in  Table  III,  a  numerical 
calculation  has  been  made  for  a  pressure  gradient  prescribed  by 


1 


1  X 
8  l 


Under  this  pressure  gradient,  we  find 


The  initial  conditions 


X 

_2  _ 

St 


c 

o 


c 

3 


are 

0.010000 

0.010006 

3.157  =  0.3157 

-1.923  =  -0.1923 

-0.3133  =  -0.3133 

1.1000 


(61) 


(62) 
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Starting  the  solution  at  4  =  0.01  rather  than  4=0  could 
influence  slightly  the  value  of  4  for  separation,  but  its  effect 
oi  the  actual  velocity  profile  near  separation  will  be  unimportant. 

The  values  of  c  ,  c  ,  c  ,  and  c  for  this  case  are  tabu- 
lated  in  Table  III  and  plotted  in  Figure  3  versus  axial  distance. 
The  separation  point  occurs  where  c3  goes  to  zero.  This  value 
of  xs/^  obtained  by  running  the  solution  on  a  fine  interval 
near  separation  is  found  to  be 

X 

-f-  =  0.9627 

This  compares  very  well  with  the  value  given  by  Leigh  of 

X 

-j-  =  0.958542  +  0.000002 

4.3  Comparison  with  Results  of  Leigh 

Leigh  has  presented  tabulated  velocity  profiles  for  small 
distances  upstream  of  separation  but  not  at  separation  itself 
since  his  numerical  methods  did  not  go  through  to  separation. 

His  positions  are  as  follows: 

Xs  X 

-j-  -  j  =  0.000142,  0.000542,  0.001542,  0.002542 

The  calculated  values  of  c  ,  c, ,  c^,  and  c  are  tabulated 
versus  X/ f  in  Table  III,  and  velocity  profiles  for  comparison 
with  those  of  Leigh  are  presented  in  Table  IV.  It  is  seen  that 
the  results  for  u  agree  with  those  of  Leigh  to  about  0.002. 

We  can  therefore  conclude  that  the  method  of  integral  relation¬ 
ships  is  yielding  excellent  approximations  to  the  true  velocity 
profiles  for  this  test  case.  We  can  also  approach  other  cases 
of  separation  with  some  degree  of  assurance  that  our  four- 
parameter  family  of  velocity  profiles  will  yield  a  close  fit. 

We  can,  of  course,  increase  the  accuracy  by  going  to  a  fifth- 
order  approximation. 
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In  comparing  our  results  with  those  of  Leigh  we  note  that  the 
coordinate  YL  of  Teigh  normal  to  the  wall  is  related  to  our 
coordinates  as  follows: 


4.4  Velocity  Profiles  at  Separation 

The  point  was  made  in  the  Introduction  that  a  one-parameter 
family  of  profiles  cannot  be  expected  to  represent  accurately  all 
possible  profiles  occurring  in  separated  and  attached  flows.  For 
instance,  one  cannot  expect  all  profiles  with  zero  wall  shear  to 
reduce  to  a  universal  function  when  considered  in  the  nondimens ional 
form  of  U/l^  versus  Y/6*.  Yet,  quartic  profiles  and  similarity 
profiles  each  yield  such  universal  functions  for  zero  shear,  and 
such  universal  profiles  would  be  expected  to  be  an  accurate  repre¬ 
sentation  of  the  separation  profile  in  all  cases  including  the 
present  calculative  example.  Since  we  have  shown  that  the  present 
calculated  profile  at  separation  is  very  accurate,  we  can  use  it 
as  a  standard  of  comparison  in  this  particular  case.  Accordingly, 
a  comparison  has  been  made  between  the  Dorodnitsyn  separation  pro¬ 
file  and  the  corresponding  quartic  and  similarity  profiles  in 
Table  V.  The  Dorodnitsyn  profile  is  accurate  to  0.002  in  U/Ui 
for  a  given  value  of  Y/5*.  The  quartic  profile  differs  from  the 
accurate  result  by  as  much  as  0.030.  The  corresponding  figure 
for  the  similarity  profile  is  0.016.  These  errors  apply,  of 
course,  only  to  the  Hartree-Leigh  case. 

The  slopes  of  the  velocity  profiles  are  compared  for  the 
present  calculative  case  and  the  quartic  case  in  Table  V(b). 
Substantial  errors  in  slope  occur  within  the  boundary  layer  for 
the  quartic  case. 
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5.  CALCULATIVE  EXAMPLE  FOR  FREE  INTERACTION 

It  is  now  of  interest  to  investigate  separated  flows  under 
the  assumption  of  free  interaction.  This  will  be  done  first  for 
the  pre-separation  region,  then  the  post-separation  region,  and 
finally  the  reattachment  region.  Some  particularly  interesting 
results  around  the  separation  point  will  be  examined.  For  the 
purposes  of  illustration  we  have  selected  a  case  for  which  data 
are  avl1  ole.  The  configuration  for  the  calculative  example  is 
showr  Figure  4. 

5.1  Initial  Conditions 

In  addition  to  the  initially  prescribed  velocity  profile,  it 
is  necessary  to  prescribe  some  initial  condition  on  pressure  for 
the  free-interact ion  equation.  The  initial  condition  should  be 
such  that  a  small  disturbance  will  initiate  the  free  interaction 
which  results  in  the  development  of  a  separated  flow.  Attempts 
were  made  to  initiate  such  an  interaction  by  perturbing  the  veloc¬ 
ity  profile,  by  specifying  a  value  of  U  ,  and  by  perturbing  p^ . 
It  was  found  that  a  small  perturbation  in  the  value  of  px  was 
a  satisfactory  way  to  induce  separation,  and  that,  as  the  per¬ 
turbation  became  small,  the  pressure  distribution  tended  to 
approach  a  limit. 

To  determine  the  change  in  Ux  accompanying  a  prescribed 
increment  in  p^  let  us  first  write  for  the  main  stream  in  the 
compressible  plane 


or 


dPl  =  -plUl  dux 


dp  du, 

1  ..  2  1 
-  =  -yM  - 

p,  1 


(64) 


We  now  utilize  the  constancy  of  the  stagnation  temperature  together 
with  the  relationship 


u ,  a 

U  - 

1  ai 


(65) 
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to  establish  the  equation 


Accordingly  we  obtain  the  desired  relationship 


dui  1  +  ^T1  Mi2  dpi 


(66) 


(67) 


giving  the  change  in  the  initial  value  of  U  to  induce  a  pre¬ 
scribed  change  dp1/p1 . 

In  addition  we  must  specify  a  value  of  0q  for  the  free- 
interaction  equation.  The  f ree-interact ion  equation  can  be  written 


tan  0 

o 


d5* 


(68) 


It  is  convenient  to  determine  this  quantity  in  the  Stewartson 
plane  where  the  displacement  thickness  5  *  and  the  momentum 

D 

thickness  6  **  are  given  for  no  pressure  gradient  by 

O 


(69) 


where  k  and  k3  are  taken  from  Reference  7.  From  Appendix  II, 
we  can  show  that,  in  general, 


o  o 


(70) 


and  in  this  case 
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b* 


(1 


+  m0)Ai 


+  m  >  „ 
o  2 


and  x  =  X. 

We  thus  have  for  the  initial  value  of  0 

o 


(71) 


(72) 


A  series  of  runs  were  made  to  see  how  the  interaction  pres¬ 
sure  distribution  behaved  as  the  magnitude  of  Apx  applied  to 
the  boundary  layer  at  xQ/i  was  reduced.  The  following  set  of 
initial  values  was  used: 


p  /p 

1  o 

1.001039 

1.01 

1.02012 

1.039251 

R 

o 

48000 

48000 

48000 

48000 

M 

O 

2.7 

2.7 

2.7 

2.7 

T 

1.4 

1.4 

1.4 

1.4 

V7 

1 

1 

1 

1 

«o 

1 

1 

1 

1 

c 

O 

3.157 

3.157 

3.157 

3.157 

c 

1 

-1.923 

-1.923 

-1.923 

-1.923 

c 

2 

-0.3133 

-0.3133 

-0.3133 

-0.3133 

c 

3 

1.1000 

1.1000 

1.1000 

1.1000 

V°o 

0.9997500 

0.9976041 

0.9952046 

0.9907388 

♦o 

0.011872 

0.011872 

0.011872 

0.011872 

M, 

2.69932 

2.69353 

2.68705 

2.67499 

The  free-interact ion  pressure  distributions  accompanying  the 
various  pressure  impulses  are  shown  in  Figure  5.  The  results 
are  readily  plotted  against  (x  -  xq)/xq  where  xq  is  the  point 


27 


AFFDL-TR-65-107 


of  the  beginning  of  interaction.  The  separation  pressure  ratio 
is  indicated  by  the  points.  It  is  clear  that  the  pressure  distri¬ 
butions  are  approximately  parallel  to  each  other  by  the  amount  of 
the  initial  impulse  with  little  change  in  the  position  01  separa¬ 
tion.  Also,  the  distributions  appear  to  tend  toward  a  limit  as 
the  impulse  approaches  zero.  Apart  from  the  question  of  the 
existence  of  a  limit,  a  pressure  impulse  of  one-tenth  of  a  per¬ 
cent  is  considered  small  for  the  present  purposes  and  is  used 
henceforth  to  initiate  interaction. 

5.2  Other  Parameters  of  Calculative  Example 

A  calculative  example  has  been  made  for  the  CS.25°-1  model 
of  Reference  8  for  which  laminar  data  are  available  in  Figure  17 
of  that  report.  The  following  initial  conditions  were  used  in 
the  calculation: 


R  =  39000,  M 
o  ’  o 


2.70,  -f  -  1.000,  5o 


1.000, 


1.001 


The  value  of  i  is  taken  as  the  distance  from  the  start  of  the 
boundary  layer  to  the  onset  of  f ree-interaction  in  the  calculative 
example.  For  this  case  we  have  taken  the  reference  length  to  be 

£ 

—  =  0.238 

4-i 

C 


where  is  given  in  Figure  4.  The  f ree-interaction  was  started 

at  the  known  experimental  location  since  the  objective  in  this 
case  is  not  to  predict  the  separation  point  from  some  reattachment 
criterion  but  to  test  the  adequacy  of  the  free-inter action  theory 
to  explain  the  observed  pressure  distribution. 

The  initial  conditions  are 


x 

-f  =  1.000,  4o  =  1.000 

c  =  3  157,  c  =  -1.923,  c,  =  -0.3133,  c  =  1.1000 

O  ’  1  »  2  *3 
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U 

“  =  0.99975,  <p  =  0.01317 
U  *  o 


The  calculation  proceeded  without  any  difficulty  and  gener¬ 
ated  the  variation  of  c  ,  c  ,  c  ,  and  c  with  x  shown  in 
Figure  6.  The  variations  downstream  of  separation  will  subse¬ 
quently  be  discussed.  The  comparison  of  the  calculated  pressures 
with  the  experimental  ones  up  to  the  separation  point  shown  in 
Figure  10  is  considered  good.  It  appears  that  the  assumption  of 
free  interaction  is  a  valid  one  up  to  separation. 

5.3  Behavior  of  the  Solution  in  the  Neighborhood  of  Separation 

One  point  of  particular  interest  is  the  nature  of  the  singu¬ 
larity  at  the  separation  point,  and  the  manner  in  which  the  free- 
interaction  theory  copes  with  this  singularity.  The  solution 
which  goes  smoothly  through  the  separation  point  is  one  based  on 
Equations  (45)  to  (48)  and  the  free-ir.teraction  relationship. 
Equations  (45)  to  (48)  have  no  discontinuities  at  separation  but 
the  free-interact ion  equation  does. 

Consider  the  free-interaction  equation 


o 


1  +  m 

o 


/ d6  *  d6.*\ 

-2-  +  — i- 

\  dx  dx  / 


R  1//2  tan  <b 

o _ _ 

1  +  m 

o 


(73) 


where  the  two  terms  on  the  left-hand  side  are  given  by  Equations 
(11-10)  and  (11-15)  of  Appendix  II.  In  crossing  the  separation 
point  the  right-hand  term  in  continuous.  Equation  (11-15)  con¬ 
tributes  a  discontinuous  term  n  on  the  left-hand  side.  How¬ 
's 

ever,  an  equal  and  opposite  discontinuity  is  contributed  in  the 

c 3 co<^o  term  of  Equation  (11-10)  so  that  equality  is  maintained. 

The  singularity  in  the  c  c  g  term  arises  from  g  since 

3  0  0  O 

%  =  2  (V1  +  s  -  Vs  )  (74) 

and,  therefore, 
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c 

3 


Since  cr  is  zero  at  separation,  g  is  infinite.  If  in  the 

3  O 

neighborhood  of  the  separation  point,  we  write 


Hs  =  \  '  V 

c3  -  V*  -  es)2 


(75) 


we  readily  obtain  across  the  separation  point 

fc  c  g  'j  -  (c  c  g  ^ 

v  3  0^0  A  _e  \  3  oyo A  + 


n  =  k  =  iim 
S  5  6-0 


£,  +e 

s 


(76) 


n  =  k  =  4c  -y  k 

's  3  o  V  4 


(77) 


We  have  thus  determined  the  slope  of  the  u  =  0  line  at  separation. 

Another  relationship  is  needed  in  addition  to  those  utilized 
upstream  to  carry  out  the  integrations  downstream  of  the  separation 
point.  It  will  now  be  shown  that  this  relationship,  compatible 
with  Equation  (77),  is  the  continuity  of  the  second  derivatives  of 
the  velocity  profiles  at  the  u  =  0  line.  If  we  use  a  quadratic 
for  the  inner  velocity  profile  we  can  expand  the  velocity  profile 
about  the  u  =  0  line  as  follows: 


u 


<T|  -  V2 


(78) 


where  the  derivatives  evaluated  from  the  outer  profile  are 


(79) 


Since  u  =  0  when  T)  =  0,  we  write 
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bif 


(80) 


T‘s  = 


^  u  ♦  ^ 

c  V  2c  c 

o  \  3  0 


The  dominant  term  as  £  — ►  £  in  this  expression  is 

s 


n  =  4c  Vc7  as  £  -*•  £ 
's  o  v  3  ! 


^S  =  4CoVk4 


(81) 


The  inner  velocity  profile  is  thus  a  quadratic  just  downstream 
of  separation,  and  the  proper  condition  required  to  continue  the 
pre-separation  solution  smoothly  through  separation  is  thus  the 
continuity  of  the  second  derivatives  of  the  inner  and  outer  veloc¬ 
ity  profiles  at  the  u  =  0  line.  Only  by  taking  the  inner 
velocity  profile  as  a  quadratic  immediately  downstream  of  separa¬ 
tion  has  it  been  possible  to  make  the  numerical  integration 
proceed. 

It  is  of  interest  to  note  that  the  conditions  of  a  quadratic 
inner  velocity  profile  just  downstream  of  separation  can  be 
inferred  from  an  analytical  continuation  argument.  Consider 
several  profiles  as  shown  in  Figure  7.  Downstream  of  separation 
the  inner  profile  corresponds  to  a  quadratic  which  will  have  two 
points  for  a  given  u  and  the  values  of  du/ch]  at  these  two 
points  will  be  equal  and  opposite.  Thus,  that  part  of  the  post¬ 
separation  profile  marked  "inner  profile"  will  hrve  the  u-axis 
as  an  axis  of  symmetry.  What  this,  in  effect,  means  is  that  the 
inner  profile  could  be  represented  as  follows: 

(1  -  u)n/u  +  c 

-  ± - = - zr  (82) 

c  +  c  u  +  c  u2 

O  1  2 


£u 
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where  the  plus  sign  refers  to  the  first  and  second  quadrant  and 
the  negative  sign  to  the  third  quadrant.  The  observed  behavior 
of  the  inner  profile  can  thus  be  explained  by  an  analytical  con¬ 
tinuation  of  the  outer  profile  into  the  inner  region. 

In  Figure  8  a  plot  is  shown  of  the  computed  values  of  c  ,  cx, 
c  ,  and  c3  in  the  neighborhood  of  the  separation  point  for  the 
case  where  the  solution  is  continued  on  the  basis  of  a  quadratic 
profile  and  a  second  derivative  match.  It  is  noted  that  Vs 
is  linear  in  x  so  that  c3  is  quadratic  as  previously  assumed. 

5.4  Adequacy  of  Boundary-Layer  Theory  in  Neighborhood  of 
Separation  Point 

While  it  has  been  shown  that  the  present  boundary-layer  cal- 
culative  method  goes  smoothly  through  the  separation  singularity, 
it  has  not  been  established  that  boundary-layer  theory  is  a 
satisfactory  representation  of  the  viscous  flow  in  such  a  neigh¬ 
borhood.  In  this  connection  we  will  show  that  the  slope  of  the 
u  =  0  line  at  separation  calculated  on  the  basis  of  the  Navier- 
Stokes  equations  is  in  good  agreement  with  the  present  calculative 
results.  We  will  also  show  that  the  use  of  free  interaction 
together  with  the  boundary-layer  equations  yields  reasonable 
physical  behavior  at  the  separation  point,  whereas  the  imposition 
of  a  prescribed  pressure  distribution  as  in  the  case  of  Leigh  and 
Hartree  can  lead  to  physically  implausible  behavior  there. 

With  regard  to  the  slope  of  the  dividing  streamline  and  the 
u  -  0  line  at  separation,  Oswatitsch  gives  the  following  results 
based  on  the  Navier-Stokes  equations: 


tan  3d 


•>  <3t/$x 
<3p/dx 


tan  0 

u 


7  d-r/dx 
dp/dx 


(83) 


The  angles  0D  and  0^  are  shown  in  Figure  9.  Based  on  the 
present  analysis  it  can  be  shown  that 
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ds 

tan  p  =  ~  = 
u  dx 


0  +  ^  “o’)  ^7 


(84) 


^  V)  ^  ft  f 


-v/c”  c 
V  3  o 


(85) 


The  following  values  at  the  first  calculated  point  beyond  separa¬ 
tion  were  found  from  the  calculations: 

M,  =  2.58,  M  =  2.70,  c  =  5.567xl0"6,  c  =  3.43X10-3 

1  ’  O  ’3  *  3 


U1  _  U! 


R  =  39000,  —  =  0.95675,  rr^  =  -0.09277,  c  =  2.190,  n  =  6.39 

o  *  U  ’  U  *  o  *  ‘s 

o  o 


The  calculated  values  from  Equations  (84)  and  (85)  are 


tan  5  =  •  0.0795 

u  dx 


-2  =  0.0854 

dp/  ox 

It  thus  appears  that  the  precise  result  of  Equation  (83)  based  on 
the  Navier-Stokes  equations  are  duplicated  within  10  percent  by 
the  boundary-layer  equations  at  the  separation  point. 

The  second  point  has  to  do  with  the  rate  of  change  of  shear 
with  axial  distance  at  the  separation  point. 


T  ~ 


br] 
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In  the  present  free-interact ion  solution  where 

c3  =  k4^  ‘  V2 

we  find  a  finite  limit  as  (•  — 1 -  £  ,  as  follows: 

s 


The  wall  skin  friction  is  thus  a  linear  function  of  £,  passing 
through  the  separation  point  as  shown  in  Figure  9.  In  the  case 
treated  by  Leigh,  c^  has  a  value  of  about  0.2  at  separation  so 
that  dx/df,  is  infinite.  The  difficulties  encountered  by 
Goldstein,  Reference  9,  in  finding  a  physically  satisfactory 
mathematical  description  of  the  flow  in  the  neighborhood  of  the 
separation  point  for  Leigh's  case  thus  appears  to  be  associated 
with  the  inadequacy  oi  the  boundary-layer  equations  in  that  case 
to  yield  a  physically  realistic  flow.  The  introduction  of  free 
interaction  into  the  problem,  however,  does  cause  the  boundary- 
layer  equation  solution  to  be  physically  realistic. 

5.5  Continuation  of  Solution  Downstream  of  Separation 

We  now  consider  continuing  the  solution  downstream  to  the 
point  where  the  ramp  starts  the  recompression.  There  is  some 
question  regarding  the  best  method  for  accomplishing  this. 

Various  cases  have  been  investigated  based  on  different  equations 
and  boundary  conditions. 

In  the  first  place,  several  different  momentum  relationships 
can  be  applied  to  the  inner  layer.  We  can  try  the  inner  momentum 
relationship  represented  by  Equation  (52).  However,  a  simpler 
wall  relationship  can  be  used  instead  based  on  applying  Equation 
(1)  to  the  wall. 

£2  *  _<L  (  ^ 

bx  by  v  dy  / 


(88) 


|y=o 
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With  regard  to  the  pressure  at  the  outer  edge  of  the  boundary  layer, 
two  possibil  ties  exist.  Firstly,  we  can  continue  to  use  the  free- 
interaction  concept  in  an  attempt  to  produce  a  unified  theory  for 
the  entire  separation  region.  Secondly,  it  is  possible  to  invoke 
the  boundary  condition  of  constant  pressure  as  exemplified  by  the 
concept  of  a  plateau  pressure  region.  Another  boundary  condition 
which  can  be  invoked  is  continuity  of  the  velocity  derivatives  at 
the  u  =  0  line. 

In  order  to  investigate  the  consequences  of  various  combina¬ 
tions  of  the  foregoing  conditions,  five  different  cases  were 
studied  as  follows: 

Case  A: 

(1)  Quadratic  inner  profile. 

(2)  Wall  relationship,  Equation  (88). 

Case  B: 

(1)  Cubic  inner  profile. 

(2)  Wall  relationship. 

(3)  Inner  momentum  equation,  Equation  (52). 

Case  C: 

(1)  Cubic  inner  profile. 

(2)  Wall  relationship. 

(3)  Continuity  of  o2u/ch)2  on  u  =  0  line. 

Case  D: 

(1)  Quadratic  inner  profile. 

(2)  Continuity  of  d2u/di)2  on  u  =  0  line. 

Case  E: 

(1)  Quadratic  inner  profile. 

(2)  Continuity  of  d2u/di,2  on  u  =  0  line. 

(3)  Constancy  of  static  pressure  from  point  where  Ui  =  0 

in  Case  D. 

Some  calculated  results  for  these  cases  are  shown  in  Figures  10 
and  11.  In  Figure  10,  the  calculated  pressure  distributions  are 
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shown  and,  ir  Figure  11,  the  calculated  velocity  profiles  and  veloc¬ 
ity  gradient  profiles.  On  the  basis  of  these  calculated  results, 
certain  conclusions  are  drawn. 

Cases  A  and  B,  while  demonstrating  reasonable  pressure  dis¬ 
tribution,  show  a  singularity  in  b2u/^rja  in  Figure  11(b)  at  the 
u  =  0  line.  Furthermore,  this  singularity  tended  to  cause  dis¬ 
turbances  in  the  velocity  profiles  which  emanate  from  the  u  =  0 
line  and  propagate  both  upward  and  downward  as  the  downstream 
distance  increases.  Such  disturbances  are  a  result  of  either 
flow  instabilities  or  calculative  instabilities;  it  is  not  known 
which.  As  a  result,  these  cases  were  discarded  and  it  was  con¬ 
cluded  that  continuity  of  d2u/cb)2  on  the  u  =  0  line  is  a 
required  boundary  condition. 

It  was  decided  to  adopt  Case  C  as  the  one  for  the  automatic 
computation  program  rather  than  Case  D  or  Case  E  because  it  uses 
a  cubic  rather  than  a  quadratic  inner  profile  and  because  it  per¬ 
mits  a  unified  f ree-interact ion  theory  for  the  entire  separated 
region.  Also,  Cases  D  and  E  do  not  satisfy  a  momentum  relation¬ 
ship  for  the  inner  flow. 

It  is  noted  that  Cases  A,  B,  and  C  involve  a  momentum  rela¬ 
tionship  for  the  inner  layer  and  that  Cases  D  and  E  do  not.  It 
is  also  noted  that  the  inner  profiles  for  A,  B,  and  C  are  much 
flatter  than  those  of  D  and  E.  An  explanation  of  this  behavior 
is  to  be  sought  in  the  momentum  equation  which  states  that  for 
constant  pressure 


db**  =  £f 
dx  2 


(89) 


If  the  wall  shear  is  forward;  that  is,  if  C^ 
should  decrease  as  x  increases.  Therefore, 
boundary-layer  profile  above  the  u  =  0  line 
that  below  it  must  flatten. 


is  negative,  6** 
if  the  part  of  the 
tends  to  fill  out, 
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5.6  Asymptotic  Behavior  of  Case  C 

No  downstream  boundary  conditions  or  reattachment  criterion 
has  yet  been  used  in  the  analysis.  The  point  of  view  is  taken 
that  a  starting  point  of  interaction  is  assumed,  and  the  downstream 
conditions  are  calculated  by  the  computer  program.  If  the  calcu¬ 
lated  downstream  conditions  are  not  in  accordance  with  some 
reattachment  criterion  (to  be  specified) ,  the  starting  point  of 
interaction  is  changed  and  a  new  set  of  downstream  conditions 
calculated.  If  a  suitable  set  of  downstream  conditions  is  found, 
then  the  solution  is  a  possible  one.  However,  a  question  of  flow 
stability  can  be  raised  at  this  point.  No  stability  analyses  are 
included  in  the  present  work. 

It  is  interesting  to  see  what  type  of  downstream  boundary 
conditions  are  developed  by  Case  C.  To  investigate  this  point, 
the  solution  was  continued  far  downstream  on  the  computing  machine. 
The  resulting  pressure  distribution  is  shown  in  Figure  12,  together 
with  calculated  positions  of  dividing  streamlines  and  u  =  0  lines. 
It  is  noted  that  as  the  pressure  falls  the  reverse  flow  behind 
separation  point  S  is  flowing  against  an  adverse  pressure  gra¬ 
dient  which  reverses  the  inner  flow.  As  a  result,  a  second  divid¬ 
ing  streamline  and  second  u  =  0  line  develop.  It  is  clear  that 
the  solutions  can  be  associated  with  two  bluff  downstream  bodies, 
one  corresponding  to  the  inner  u  =  0  line  with  a  prescribed 
normal  velocity  distribution  and  the  other  corresponding  to  the 
inner  dividing  streamline  with  a  prescribed  slip  condition.  In 
the  concept  of  free  interaction,  the  pressure  distribution  up  to 
the  plateau  region  is  independent  of  the  downstream  means  of 
causing  separation.  (The  pressure  distribution  starting  at  the 
compression  surface,  of  course,  does  depend  on  the  shape  of  the 
compression  surface.)  Because  separation  usually  occurs  close 
to  the  compression  surface,  no  appreciable  falling  pressure  occurs 
in  the  plateau  region.  It  thus  appears  unnecessary  to  invoke  a 
constant  plateau  pressure  as  a  boundary  condition. 

The  foregoing  flow  field  actually  corresponds  to  two  stand¬ 
ing  vortices  in  the  separated  region  with  reversed  signs.  To 
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obtain  this  pattern  requires  going  to  a  cubic  inner  profile.  It 
is  not  clear  if  going  to  a  quar  .ic  profile  will  produce  three 
standing  vortices.  We  have  certainly  not  developed  a  convergent 
solution  for  the  inner  flow,  although  the  outer  flow  should  be 
quite  accurate  for  the  solid  boundary  condition  previously  men¬ 
tioned.  Higher  approximations  can  be  made  by  the  present  method. 
However,  if  a  second  separation  point  developes,  we  can  argue 
logically  that  the  inner  profile  should  be  represented  by  a  veloc¬ 
ity  profile  with  a  square  root  singularity  comparable  to  that  for 
the  outer  flow  based  on  the  same  arguments  advanced  for  the  square 
root  singularity  in  Equation  (36) . 

Let  us  examine  the  velocity  profile  of  the  outer  flow  to  see 
if  an  asymptotic  outer  profile  is  developed.  It  turns  out  that 
the  profile  outside  the  outer  u  =  0  line  approaches  that  given 
by  Chapman  in  Reference  10.  Also,  the  velocity  profile  can  be 
non-dimensionalized  with  respect  to  x  in  the  same  manner  as  the 
Chapman  solution.  In  this  form  the  velocity  profile  U/l^  is  a 
function  only  of  the  parameter  £  defined  by  Chapman  as  follows: 

y 

if  pudy 

- - -  (90) 
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The  correspondence  between  the  foregoing  notation  of  Chapman  and 
the  notation  of  the  present  report  is 


u 


=  u 


=  V 


T  =  T 

00  l 


qi  I 


In  the  present  analysis  a  parameter  f  is  evalua ted  as  follows: 


(92) 
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It  is  then  easy  to  show  that 


so  that  a  direct  comparison  with  the  results  of  Chapman  can  now 
be  made . 

A  comparison  is  made  between  the  velocity  profiles  of  Case  C 
and  that  of  Chapman  (Ref.  10)  in  Figure  13.  The  two  points  where 
ip  =  0  and  the  two  points  where  u  =  0  are  seen  in  the  velocity 
profile  for  Case  C.  The  velocity  profile  above  the  u  =  0  line 
originating  at  S'  shows  a  close  approximation  to  that  of  Chapman 
in  an  asymptotic  sense. 

5.7  Continuation  of  Solution  to  Reattachment  Point 

The  present  analysis  allows  for  introduction  of  a  compression 
ramp  of  fairly  arbitrary  shape  at  some  point  downstream  of  separa¬ 
tion.  Case  C  has  been  made  to  go  through  to  the  reattachment 
point  by  observing  that  the  velocity  profile  approaching  the 
reattachment  point  is  quadratic  just  as  is  the  profile  coming  out 
of  the  separation  point.  As  the  solution  for  Case  C  proceeds 
downstream  of  the  beginning  of  the  ramp,  the  cubic  velocity  pro¬ 
file  becomes  quadratic.  At  that  point,  a  quadratic  profile  is 
assumed  and  continued  up  to  reattachment.  In  this  way  the  solu¬ 
tion  was  made  to  go  smoothly  into  the  reattachment  point. 

In  Figure  14,  the  calculated  flow  characteristics  are  shown. 
The  value  of  c3  and  c3  are  both  zero  at  the  reattachment  point 
just  as  at  the  separation  point.  At  a  point  about  half  way  up 
the  compression  surface  toward  the  reattachment  point,  the  cubic 
profile  turns  into  a  quadratic  profile.  At  this  point  the  pro¬ 
file  is  taken  as  quadratic  and  is  so  continued  into  the  reattach¬ 
ment  point.  The  dividing  streamline  is  shown  together  with  the 
u  =  0  line  which  is  about  two-thirds  of  the  distance  between 
the  wall  and  the  dividing  streamline. 
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Comparison  of  the  calculated  and  experimental  pressure  dis¬ 
tributions  shows  a  long  delay  downstream  of  the  beginning  of  the 
ramp  before  the  pressure  starts  its  rise.  The  significance  of 
this  fact  can  be  seen  by  comparison  with  the  inviscid  pressure 
distribution.  The  calculative  and  experimental  pressure  distribu¬ 
tions  are  in  fair  agreement.  The  reversal  in  curvature  of  the 
calculated  pressure  distribution  approaching  reattachment  is 
exaggerated  by  the  large  scale  and  is  not  nearly  so  pronounced  in 
subsequent  calcu lat xcns  to  be  shown  for  wedges.  It  thus  appears 
to  be  dae  to  the  increasing  slope  of  the  ramp.  As  a  matter  of 
interest,  the  calculated  total  pressure  on  the  dividing  streamline 
is  also  shown. 

In  Figure  12  it  was  seen  that  without  a  compression  ramp, 
the  reversed  flow  separated  at  a  point  S'  located  at  x/ £  of 
3.9.  Since  the  compression  ramp  is  introduced  at  x/i  of  about 
2.1,  we  will  still  have  separation  at  S'  with  tne  ramp.  The 
calculated  separation  region  and  velocity  pi.  Les  for  Case  C 
with  the  ramp  are  illustrated  in  Figure  15.  It  is  seen  that  the 
ramp  causes  the  inner  separated  flow  to  reattach  very  rapidly. 

In  fact,  the  small  amount  of  inner  separated  flow  and  its  very 
low  velocity  makes  it  of  no  physical  importance  insofar  as  the 
overall  flow  is  concerned. 
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6.  CALCULATIVE  PROGRAM 

6.1  General  Aspects  of  Calculative  Program 

It  has  been  the  primary  objective  of  this  study  to  develop 
an  automatic  computational  program  for  calculating  laminar  sepa¬ 
rated  flows.  The  program  which  has  been  developed  is  described 
completely  in  an  operating  manual.  While  it  is  not  intended  to 
describe  in  any  detail  the  calculative  program  herein,  neverthe¬ 
less,  certain  general  features  of  the  program  are  of  interest 
and  should  aid  in  utilizing  the  program  to  maximum  advantage. 

The  calculative  technique  is  basically  an  iterative  one .  The 
beginning  of  interaction  between  the  boundary  layer  and  the  outer 
flow  is  assumed  to  start  at  a  certain  position,  and  the  resulting 
separated  flow  is  calculated  up  to  the  reattachment  point.  If 
the  calculated  reattachment  point  does  not  fulfill  the  particular 
reattachment  criterion  which  may  be  invoked,  one  then  changes  the 
beginning  of  interaction  and  computes  a  new  reattachment  point. 

In  this  sense  the  method  is  an  iterative  one.  No  particular 
reattachment  criterion  has  been  built  into  the  method.  No  guaran¬ 
tee  of  the  existence  of  a  solution  is  given.  Even  if  a  solution 
is  found  by  calculation,  no  guarantee  of  its  stability  is  made. 

A  discussion  of  the  reattachment  criterion  will  subsequently  be 
made . 

The  method  can  be  applied  to  a  compression  or  expansion  ramp 
of  fairly  general  shape,  but  the  particular  program  here  is  written 
for  a  flat  plate  with  a  flat  wedge  faired  together  by  a  circular 
arc.  The  radius  of  this  arc  can  be  varied.  Within  reasonable 
limits,  the  radius  size  does  not  influence  the  calculated  results. 

The  velocity  profile  external  tD  the  u  =  0  line  has  been 
taken  as  a  four-parameter  profile,  and  the  inner  profile  has  been 
taken  as  a  cubic  except  near  separation  and  reattachment.  On  the 
basis  of  comparison  of  calculated  results  by  the  present  method 
with  "exact  solutions"  from  other  methods  and  on  the  basis  of  com¬ 
parison  between  calculated  results  and  experimental  results,  it 
is  believed  that  the  above  profiles  will  give  sufficiently  accurate 
approximations  to  the  flow  field  for  most  purposes.  Possibilities 
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for  higher  approximations  are  inherent  in  the  method,  but  none  has 
been  carried  out.  This  should  be  done. 

The  program  is  based  on  the  concept  of  free  interaction  between 
the  boundary  layer  and  the  external  flow  in  that  the  pressure  dis¬ 
tribution  acting  on  the  body  in  viscous  flow  is  assumed  to  be  the 
same  as  the  inviscid  pressure  distribution  for  the  body  thickened 
locally  by  the  amount  of  the  boundary-layer  displacement  thickness. 
The  Prandtl-Meyer  relationship  is  used  to  calculate  the  pressure. 

The  calculative  method  will  be  accurate  to  the  extent  that  "free 
interaction"  accurately  describes  the  fluid  mechanics  of  the  separa¬ 
tion  and  reattachment  processes.  The  calculation  is  started  by 
assuming  that  the  flow  up  to  some  point  termed  the  beginning  of 
interaction  is  of  the  Blasius  type,  and  that  free  interaction 
starts  discontinuous ly  at  that  point.  The  precise  formulation  of 
the  initial  conditions  will  be  subsequently  discussed.  The  intro¬ 
duction  of  free  interaction  together  with  a  small  pressure  impulse 
applied  to  the  boundary  layer  is  sufficient  to  cause  the  computing 
scheme  to  generate  the  separated  flow  solution. 

From  the  mathematical  point  of  view  the  computational  program 
is  simply  one  that  solves  a  simultaneous  set  of  ordinary  differen¬ 
tial  equations  step  by  step  starting  at  the  beginning  of  interac¬ 
tion  and  continuing  downstream  to  the  reattachment  point.  The 
set  of  equations  is  obtained  from  the  method  of  integral  rela¬ 
tions  applied  to  the  laminar-boundary-layer  equations  and  from 
the  free-interaction  relationship.  The  unknowns  in  the  equations 
are  the  parameters  which  specify  the  velocity  profiles  and  the 
pressure  distribution.  For  the  region  in  front  of  the  separation 
point  there  are  four  velocity  profile  parameters  which,  taken 
together  with  the  pressure,  require  five  simultaneous  differential 
equations  for  their  determination.  Downstream  of  the  separation 
point  a  cubic  profile  is  utilized  for  the  inner  flow  with  three 
parameters  determining  the  inner  velocity  profile,  Equation  (38). 
However,  the  use  of  the  boundary  condition  of  the  continuity  of 
two  velocity  derivatives  at  the  boundary  between  the  inner  and 
the  outer  flows  results  in  the  addition  of  only  one  more  differ¬ 
ential  equation  downstream  of  the  separation  point. 
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The  program  can  be  utilized  for  arbitrary  values  at  the 
beginning  of  interaction  of  the  Reynolds  number,  Mach  number,  and 
ratio  of  specific  heats.  By  eliminating  the  equation  of  free 
interaction  and  specifying  an  external  pressure  distribution,  the 
method  can  be  used  to  calculate  with  excellent  accuracy  the  change 
in  shape  of  the  laminar  velocity  profile  in  the  face  of  an  arbi¬ 
trary  pressure  gradient.  It  is  possible  to  increase  the  accuracy 
of  the  program  by  increasing  the  order  of  approximation  used  to 
represent  the  velocity  profiles.  (Until  the  stability  of  the 
calculative  process  has  been  demonstrated,  no  positive  claims  can 
be  made.)  It  is  also  capable  of  extension  to  much  broader  appli¬ 
cation  as  discussed  in  "Recommendations  for  Future  Work." 

6.2  Assumptions  of  Calculative  Program 

A  number  of  assumptions  have  been  made  in  the  analysis  as 
listed  in  Section  2.2.  In  addition,  several  additional  assump¬ 
tions  or  conditions  have  been  specified  for  the  calculative  pro¬ 
gram. 

(1)  A  four-parameter  velocity  profile  has  been  used  for  the 
outer  flow. 

(2)  A  cubic  profile  has  been  used  for  the  inner  flow. 

(3)  The  configuration  is  a  flat  plate  with  a  wedge  com¬ 
pression  surface  joined  by  an  arc  of  constant  radius. 

(4)  The  initial  conditions  of  the  boundary  layer  correspond 
to  those  for  the  Blasius  case. 

6.3  Initial  Conditions 

Some  discussion  of  initial  conditions  is  required  for  cases 
of  very  low  Reynolds  numbers  or  very  large  Mach  numbers .  For 
these  cases  it  is  found  that  the  streamwise  slope  of  the  displace¬ 
ment  thickness  at  the  beginning  of  interaction  is  not  small,  being 
of  the  order  of  10°  in  some  cases.  In  the  work  contained  through 
the  first  five  sections,  the  assumption  was  made  in  calculating 
the  initial  conditions  that  the  velocity  ui  at  the  edge  of  the 
boundary  layer  is  parallel  to  the  plate.  However,  it  was  found 
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necessary  to  take  account  of  the  difference  in  direction  between 
the  velocity  at  the  edge  of  the  boundary  layer  and  the  direction 
parallel  to  the  plate  for  the  above  cases  if  the  calculative 
solution  is  to  start  smoothly. 

Consider  the  following  sketch  which  illustrates  the  initial 
condition  problem: 


At  some  point  xq  we  assume  that  the  static  pressure  jumps  dis¬ 
continuous  ly  from  pQ  to  px.  This  increment  is  small,  being 
only  one-tenth  of  1  percent  of  pQ.  Up  to  xq  one  assumes  the 
pressure  was  constant  and  that  the  boundary-layer  profile  just 
before  xq  corresponds  to  that  for  uniform  pressure.  In  some 
cases,  interaction  between  the  boundary  layer  and  the  outer  flow 
will  produce  large  induced  pressure  gradients  before  xq.  We 
have  neglected  such  pressure  gradients  or  alternately  consider 
that  the  plate  was  cambered  to  counteract  them  and  to  produce  a 
uniform  pressure.  At  xq  the  slope  of  the  6*  curve  changes 
by  a  very  small  amount  from  0q  to  01  in  accordance  with  the 
Prandtl-Meyer  relationship.  The  Mach  numbers  Mq  and  Mx  are 
defined  on  the  basis  of  the  velocity  parallel  to  the  wall. 
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Applying  the  Prandtl-Meyer  relationship  to  the  slope  of  the 
streamline  at  5*  with  the  outer  flow  Mach  number  defined  parallel 
to  the  plate,  we  find 


( 


(94) 


We  wish  to  determine  the  values  of  Mx  and  0i  as  initial 
conditions  for  the  calculation  program.  In  the  preceding 
equation  we  have  used  the  following  initial  values  to  calcu 
late  M  : 


1.001 


(95) 


(96) 


The  second  relationship  is  based  on  Equation  (72)  and  neglects  the 
small  difference  between  0  and  0,  in  the  calculation  of  M, . 

O  x  ± 

However,  in  determining  the  initial  value  of  0  itself,  we 
utilize  the  relationship  between  Mach  number  and  turning  angle. 


(97) 
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The  constant  C  is  evaluated  from  the  conditions 

M  =  M  0=0 

o  o 

so  that 

tan'1  V(^)2 

(98) 

We  then  compute  0  from  Equation  (97)  using  C  evaluated  from 
Equation  (98)  with  M/cos  0  equal  to  M1/cos  0  .  Iteration  of 
Equations  (94)  to  (97)  to  establish  more  precise  values  of  Mi 
and  0i  is  usually  not  required  because  of  the  small  difference 
between  0  and  0  . 

In  the  calculations  discussed  through  Section  5,  the  cos  0 
term  was  included  in  the  pressure  distribution  calculations  but 
was  not  included  in  the  initial  condition  calculations.  It  has, 
however,  been  included  in  both  calculations  henceforth. 

6.4  Calculative  Accuracy 

In  order  to  assess  to  what  extent  any  differences  between 
experiment  and  theory  can  be  ascribed  to  theory  rather  than  inac¬ 
curacies  in  numerical  methods,  it  is  important  to  investigate 
the  calculative  accuracy  of  the  method.  In  this  connection,  we 
will  consider  the  question  of  interval  size,  corner  radius,  and 
order  of  approximation  of  the  velocity  profiles. 

First,  a  point  concerning  stability  of  the  calculative  method 
should  be  mentioned.  If  a  laminar  boundary  layer  on  a  flat  plate 
is  given  a  distorted  profile  as  an  initial  condition  and  is  per¬ 
mitted  to  proceed  along  the  plate  at  uniform  pressure,  it  will 
develop  a  Blasius  profile  asymptotically  to  any  desired  numerical 
accuracy.  If  we  were  to  assume  a  Blasius  profile  at  some  point 
and  were  to  integrate  upstream  to  find  the  initial  conditions,  we 
would  have  a  numerically  intractable  problem  since  any  number  of 
initial  conditions  could  result  in  a  final  Blasius  profile  within 
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prescribed  numerical  accuracy.  Integrating  the  boundary-layer 
equations  downstream  of  separation  into  the  reversed  inner  flow 
has  some  of  the  aspects  of  the  foregoing  problem.  In  this  case, 
however,  it  is  complicated  by  the  existence  of  inner  and  outer 
flows.  It  was  found  that  if  the  inner  profile  was  represented  as 
a  cubic,  the  integration  immediately  downstream  of  separation 
became  unstable.  In  the  program,  the  inner  profile  is  started  as 
a  quadratic  and  changed  to  cubic  when  numerical  differences  of 
the  order  10  4  between  the  magnitude  of  the  slopes  at  the  wall 
and  the  u  =  0  line  is  found  in  a  side  calculation.  This  pro¬ 
cedure  eliminates  the  instability  completely.  It  is  not  known 
what  caused  the  instability  for  the  cubic  case.  The  other  insta¬ 
bility  encountered  in  the  calculations  occurred  as  a  result  of  a 
singularity  in  the  second  derivative  of  the  inner  and  outer  veloc¬ 
ity  profiles  at  their  juncture. 

The  interval  size  is  not  constant  in  the  integration  scheme 
but  is  variable  to  maintain  certain  prescribed  accuracies.  At 
the  start  of  the  numerical  solution  a  Runge-Kutta  integration 
scheme  is  used  for  the  first  four  steps.  Thereafter  a  fourth- 
order  Adams  predictor-corrector  method  is  used.  With  the  value 
of  the  variables  known  for  steps  j-3,  j  -  2,  j  -  1,  and  j, 
the  value  for  step  j  +  1  is  predicted  by  the  following  formula: 


B-?!  =  B.  +  -^7  ( 55B.  -  596.  +  37B .  _  -  9B .  ^  (99) 

3+i  3  24  V  3  l-i  I”2  1-3 ) 


Using  the  predicted  value  of 
is  calculated  as  follows: 


•  (  n ) 

,  a  corrected  value  of 
l+i* 


(c 

j 


b!^  =  B.  +  (9b!pJ  +  19B.  -  5B.  +  B.  S) 

3+1  3  24  V  ]H  3  I"1  l-2y 


(100) 


The  truncation  error  Ax  is  estimated  by  the  formula 


Ax  =  -  B (c)  -  B(p)  ^ 

270  l  3+1  j+i/ 


(101) 


The  absolute  error  estimate  is 
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(c) 


14.2 


(102) 


and  the  relative  error  estimate  is 


JM 

R(C) 

Vi 


(103) 


For  the  calculations  of  this  report,  the  relative  error  has 
been  taken  as  10  5  and  the  absolute  error  as  10  7 .  If  the  rela¬ 
tive  error  exceeds  10  5  but  the  absolute  error  is  less  than  10  7, 
the  integration  proceeds.  If  both  tests  fail,  the  increment  is 
halved,  and  the  calculation  started  again  using  the  Runge-Kutta. 

If  the  relative  error  is  less  than  1/200  of  the  relative  error 
bound,  the  interval  is  doubled.  Also,  if  the  absolute  error  is 
less  than  1/200  of  the  absolute  error  bound,  the  increment  is 
doubled  even  though  the  relative  error  test  fails.  These  criteria 
are  established  to  avoid  oscillating  between  two  interval  sizes. 
These  tests  are  built  into  the  IBM  Fortran  II  subroutine  used  to 
solve  the  system  of  differential  equations.  A  discussion  of  the 
method  and  the  Runge-Kutta  starting  technique  is  presented  ir 
F.  D.  Hildebrand's  "Introduction  to  Numerical  Analysis,"  McGraw-Hill, 
New  York,  1956,  Chapter  6. 

During  the  calculations  run  for  comparison  with  the  results 
of  Hartree  and  Leigh,  the  allowable  relative  error  was  decreased 
from  10  5  to  10  4 .  The  calculated  values  of  x  for  a  given  value 
of  c3  did  not  change  by  more  than  one  in  10  4.  This  was  also 
the  accuracy  of  prediction  of  the  position  of  the  separation  point. 
Accordingly,  a  value  of  relative  error  of  10  5  was  considered 
sufficiently  accurate  for  the  calculations. 

In  order  not  to  introduce  any  discontinuities  which  might 
generate  calculative  instabilities,  it  was  decided  to  fair  the 
plate  into  the  ramp  with  a  circular  arc.  For  this  reason,  calcu¬ 
lations  were  made  for  the  configuration  of  Figure  16  with  fairings 
of  several  radii.  The  calculated  results  for  the  pressure  distri¬ 
bution  are  shown  in  Figure  17.  It  is  seen  that  fairing  the  corner 
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causes  modification  of  the  pressure  distribution  only  over  its 
length  and  that  the  final  pressure  distribution  on  the  ramp  does 
not  depend  on  the  size  of  the  radius.  Also,  as  the  size  of  the 
radius  tends  to  zero,  the  calculated  pressure  distributions  tend 
to  a  limit.  It  thus  appears  that  fairing  the  corner  causes  no 
inaccuracies  in  the  calculated  pressure  distributions.  Some 
1  "»cal  modifications  of  the  velocity  profiles  may  be  anticipated 
in  the  region  of  the  fairing.  In  the  program,  the  actual  radius 
of  the  fairing  is  left  arbitrary. 

It  has  been  mentioned  that  the  present  method  contains  within 
its  framework  the  possibility  of  obtaining  higher-order  approxima¬ 
tions  to  the  calculated  flow.  Such  approximations  are  obtained  by 
increasing  the  number  of  parameters  in  the  velocity  profiles.  With 
regard  to  the  outer  profile,  a  four-parameter  formulation  has  been 
used  throughout  the  paper  and  a  higher-order  approximation  has  not 
been  investigated.  Based  on  the  excellent  agreement  of  the  four- 
parameter  approximation  with  the  results  of  Leigh,  it  was  felt 
that  a  fourth  approximation  for  the  outer  flow  is  enough  for  most 
purposes.  Furthermore,  we  must  consider  the  question  of  increas¬ 
ing  the  degree  of  approximation  of  the  outer  profile  simultaneously 
with  that  for  the  inner  profile. 

It  is  felt,  insofar  as  the  pressure  distribution  is  concerned, 
the  number  of  arbitrary  parameters  in  the  inner  profile  can  be 
less  than  that  for  the  outer  profile.  In  the  case  of  the  cubic 
profile,  we  have  matched  the  velocity  and  its  first  two  derivatives 
at  the  boundary  between  the  inner  and  outer  flow,  leaving  only 
one  free  parameter  to  satisfy  an  inner  momentum  relationship 
(Case  C) .  For  Case  D  with  a  quadratic  inner  profile,  continuity 
of  the  velocity  and  its  first  two  derivatives  fully  determines 
the  profile,  and  momentum  is  not  satisfied  for  the  inner  layer 
(although  continuity  is) . 

A  comparison  between  Cases  C  and  D  and  the  data  for  the  par¬ 
ticular  case  is  shown  in  Figure  18.  It  is  clear  from  a  comparison 
of  the  calculated  results  for  Cases  C  and  D  that  for  accurate 
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results  we  must  use  a  cubic  inner  profile.  For  the  cubic  profile, 
the  agreement  between  experiment  and  theory  is  considered  fairly 
good.  How  much,  it  would  be  further  improved  by  going  to  a  quartic 
inner  profile  is  a  question  recommended  for  future  study. 
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7.  SOME  CALCULATIVE  EXAMPLES 
7.1  Systematic  Calculations 

In  order  to  check  out  the  calculative  program  over  a  wide 
range  of  Reynolds  numbers,  Mach  numbers,  and  ratios  of  specific 
heat,  a  systematic  series  of  computation  runs  were  made  by  the 
IBM  7094  Computer  at  RTD,  W-PAF3,  under  the  supervision  of 
Mr.  Eugene  Fleeman,  Project  Monitor.  The  following  tabulated  sets 
are  presented  to  show  the  influence  on  the  above  parameters  for  a 
flat  plate  with  a  25°  compression  wedge: 
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A  value  of  00  for  x_  indicates  that  reattachment  did  not  take 
place. 

The  systematic  effect  of  Reynolds  number,  Mach  number,  and 
ratio  of  specific  heats  is  shown  in  Figure  19.  In  this  and  sub¬ 
sequent  figures,  the  reference  length  &  is  taken  equal  to  xq. 

In  fact,  this  equality  is  built  into  the  program  so  that  inter¬ 
action  begins  when  x/i  is  unity. 

The  results  for  the  effect  of  Reynolds  number  show  a  decrease 
in  the  distance  from  the  onset  of  interaction  to  separation  as  the 
Reynolds  number  increases,  as  well  as  a  decrease  in  the  value  of 
the  plateau  pressure.  The  results  cannot  be  compared  directly 
with  the  earlier  results  of  this  study  in  Figure  2(a)  of 
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Reference  1  because  the  beginning  of  interaction  has  been  defined 
differently  in  each  case.  However,  the  plateau  pressure  levels 
should  be  comparable  and  do  show  good  agreement. 

With  regard  to  Mach  number,  there  is  an  increase  in  plateau 
pressure  and  a  decrease  in  the  distance  between  the  beginning  of 
interaction  and  the  location  of  separation  as  the  Mach  number 
increases.  The  plateau  pressures  for  Mach  numbers  of  2,  4,  and 
6  are  in  good  accord  with  the  earlier  calculated  values  in  Fig¬ 
ure  2(b)  of  Reference  1.  As  the  ratio  of  specific  heats  increases, 
Figure  19(c)  shows  an  increase  in  the  distance  to  the  separation 
point  but  no  significant  change  in  the  plateau  pressure. 

7.2  Comparison  Between  Theory  and  Experiment 

The  calculative  method  has  been  applied  to  several  cases 
for  which  data  on  laminar  separation  and  reattachment  are  avail¬ 
able  in  order  to  assess  its  accuracy.  The  Ames  Research  Center 
of  NASA  has  furnished  several  sets  of  such  data  for  a  flat  plate 
with  a  10°  wedge  placed  at  various  angles  in  the  wind  tunnel. 

These  data  and  the  calculated  pressure  distributions  are  compared 
in  Figures  20  and  21. 

In  obtaining  the  calculated  pressure  distributions,  the 
beginning  of  interaction  was  varied  until  the  position  of  the 
calculated  separation  point  coincided  as  closely  as  could  be 
estimated  with  the  experimental  position  of  the  separation  point. 
There  was  relatively  little  change  in  separation  point  as  a 
result  of  changes  in  plate  angle  of  attack.  Comparison  is  shown 
for  the  upper  surface  for  angles  of  attack  of  6°,  15°,  and  16.7°. 

It  is  felt  that  the  agreement  between  the  experimental  and  cal¬ 
culated  results  is  good.  It  is  of  interest  to  speculate  whether 
the  agreement  would  be  further  improved  by  going  to  a  higher 
approximation  in  the  calculative  method.  Also,  it  is  of  interest 
to  inquire  what  part  of  the  difference  between  experiment  and 
theory  can  be  ascribed  to  the  assumption  of  free  interaction. 

These  questions  are  a  fitting  subject  for  future  investigation. 
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8.  FLOW  IN  NEIGHBORHOOD  OF  REATTACHMENT  POINT 

Let  us  now  use  the  calculative  program  to  study  the  nature 
of  the  flow  in  the  neighborhood  of  the  reattachm^nt  point  and  to 
study  relationships  between  the  locations  of  separation  and 
reattachment.  First,  a  series  of  computations  were  made  for  the 
flat  plate  with  a  10°  wedge  systematically  varying  the  position 
of  the  beginning  of  interaction.  The  calculated  pressure  distri¬ 
butions,  separation  points,  and  reattachment  points  are  shown  in 
Figure  21  where  they  are  compared  with  the  measured  pressure  dis¬ 
tribution.  Two  points  of  interest  are  immediately  evident.  For 
all  the  cases  calculated,  there  are  no  substantial  variations 
among  the  calculated  wedge  pressure  distributions;  that  is,  the 
position  of  reattachment  does  not  seem  to  influence  the  wedge 
pressure  distribution  substantially.  Secondly,  small  changes  in 
the  position  of  the  beginning  of  interaction  are  sufficient  to 
cause  large  movements  in  the  reattachment  point  and,  for  the  most 
upstream  position,  reattachment  did  not  occur  at  all. 

It  is  of  interest  to  examine  the  slope  of  the  dividing  stream¬ 
line  approaching  reattachment  and  to  see  how  this  is  related  to 
the  calculated  Oswatitsch  reattachment  angle,  Equation  (83).  For 
this  purpose,  the  calculated  solution  from  the  previous  figure 
corresponding  to  i  -  0.0550  was  used  with  the  results  shown  in 
Figure  22 .  It  is  noted  that  of  the  dividing  streamline 

decreases  approaching  the  reattachment  point.  The  Oswatitsch 
reattachment  angle  calculated  from  Equation  (83)  using  the  numeri¬ 
cal  values  for  the  present  case  increases  and  then  decreases 
approaching  reattachraent .  The  values  at  the  reattachment  posi¬ 
tion  are  considered  to  be  in  good  accord  considering  one  theory 
is  based  on  the  boundary-layer  equations  and  the  other  on  the 
Navier-Stokes  equations.  For  position  of  reattachment  lower  on 
the  wedge,  the  agreement  is  not  as  good. 

Examination  of  the  calculated  Oswatitsch  reattachment  angle 
as  it  varies  with  the  position  of  the  beginning  of  interaction  l 
shows  that  the  reattachment  point  has  an  upper  limiting  position 
on  the  wedge.  To  show  this  result.  Figure  23  has  been  constructed. 
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The  variation  of  Oswatitsch  reattachment  angle  for  the  previous 
figure  is  shown  together  with  those  for  the  other  positions  of 
the  beginning  of  interaction.  It  is  seen  that  if  the  beginning 
of  interaction  starts  at  i  =  0.0525,  reattachment  cannot  occur 
for  a  greater  value  of  x  from  0.325  feet.  At  this  point  the 
dividing  streamline  is  tangent  to  the  wedge. 

On  the  basis  of  the  present  calculations,  it  is  seen  that 
the  reattachment  angle  decreases  as  the  reattachment  point  moves 
up  the  wedge.  An  extrapolation  of  the  calculated  results  to 
of  zero  yields  a  value  of  x  =  0.34.  The  present  calculation 
thus  yields  the  same  result  as  the  Oswatitsch  theory  in  that  the 
flow  cannot  reattach  to  the  wedge  above  a  certain  point.  This 
condition  corresponds  to  the  most  upstream  allowable  position 
of  the  beginning  of  interaction  and  to  a  dividing  streamline 
tangent  to  the  wedge  at  reattachment . 

It  has  been  stated  that  the  prerent  calculative  scheme  is  an 
iterative  one  in  that  the  separation  point  can  be  moved  until 
some  reattachment  criterion  can  be  fulfilled.  Insofar  as  the 
present  calculative  method  is  concerned,  reattached  flows  are 
possible  for  a  certain  range  of  reattachment  on  the  wedge.  Which 
one  of  these  flows  will  occur  in  practice  depends  on  downstream 
boundary  conditions,  and  in  this  sense  the  present  calculative 
scheme  is  an  elliptical  one.  One  speculates  that  if  the  length 
of  the  wedge  is  less  than  the  length  over  which  reattachment  will 
occur  that  reattachment  might  occur  at  the  end  of  the  wedge. 

This  speculation  is  based  on  the  physical  argument  that  the  end 
of  the  wedge  is  the  only  distinguished  point.  If  the  wedge  is 
longer  than  the  critical  value,  one  speculates  that  reattachment 
might  occur  at  the  critical  length  from  the  corner  of  the  wedge 
since  this  position  is  a  distinguished  one.  These  speculations 
require  a  stability  analysis  of  the  flow  to  establish  their 
validity.  In  addition  to  a  stability  analysis,  careful  experi¬ 
mental  measurements  are  required. 

As  a  point  of  interest,  the  variations  in  stagnation  pres¬ 
sure  and  Mach  number  along  the  dividing  streamline  for  the  several 
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cases  are  shown  in  Figure  24.  At  reattachment,  the  stagnation 
pressure  on  the  dividing  streamline  is  equal  to  the  static  pres¬ 
sure  at  this  point  and  these  points  are  shown  as  solid  points  in 
Figure  24(a).  It  is  noted  that  for  the  case  that  did  not  reattach, 
the  stagnation  pressure  continues  to  rise.  The  Mach  number  shows 
an  almost  linear  decrease  with  distance  as  it  approaches  reattach¬ 
ment  . 
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9.  CONCLUSIONS 

The  present  paper  presents  a  method  for  calculating  the 
laminar  compressible  separated  flow  past  a  flat  plate  with  a 
ramp  of  arbitrary  shape  and  small  slopes  from  an  assumed  position 
for  the  beginning  of  interaction  to  the  reattachment  point.  The 
method  is  based  on  the  assumption  of  free  interaction  throughout 
the  separated  flow  region. 

(1)  The  method  has  been  successfully  programed  for  a  flat 
plate  with  a  flat  ramp  connected  by  means  of  a  circular  arc. 

(2)  For  the  several  cases  investigated,  reasonably  good 
agreement  was  obtained  between  the  experimental  and  theoretical 
pressure  distributions  up  to  the  reattachment  point. 

(3)  In  the  neighborhood  of  the  separation  point  induced  by 
a  linear  decreasing  velocity,  the  present  method  yielded  numeri¬ 
cal  results  in  close  agreement  with  the  precise  results  of  Hartree 
and  Leigh  based  on  finite  difference  equations. 

(4)  For  an  arbitrarily  imposed  pressure  distribution,  there 
appears  to  be  a  singularity  in  the  rate  of  change  of  wall  skin 
friction  with  distance  along  the  wall.  The  assumption  of  free 
interaction  between  the  outer  flow  and  the  boundary  layer  removes 
this  singularity.  On  this  basis  the  present  solution  was  made  to 
go  smoothly  through  the  separation  point. 

(5)  The  calculated  slope  of  the  dividing  streamline  at 
separation  is  in  good  accord  with  the  value  due  to  Oswatitsch  as 
calculated  from  the  Navier-Stokes  equations. 

(6)  The  analysis  indicates  that  the  velocity  profile  just 
downstream  of  separation  and  just  upstream  of  reattachment  is 
basically  quadratic. 

(7)  The  calculative  method  yielded  reversal  of  the  inner 
flow  in  certain  cases. 

(8)  It  was  found  for  a  specific  case  that  reattached  flow 
was  possible  only  over  a  region  of  the  wedge  within  a  critical 
distance  from  the  lower  corner.  At  this  critical  distance,  the 
dividing  streamline  comes  in  tangent  to  the  wedge  and  the  separa¬ 
tion  point  is  at  its  most  forward  position. 
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(9)  The  method  contains  the  possibility  of  high  orders  of 
approximation  and  of  extension  to  other  flow  problems  as  described 
in  the  next  section. 

(10)  The  calculative  program  has  been  applied  to  a  wide  range 
of  initial  Mach  numbers,  Reynolds  numbers,  and  ratio  of  specific 
heats  to  show  the  systematic  effects  of  these  parameters  on  the 
separated  flow. 

(11)  The  present  method  can  be  applied  to  the  calculation 
of  the  laminar -boundary-layer  characteristics  in  the  face  of  a 
proscribed  pressure  distribution  for  any  Mach  number. 

(12)  It  is  believed  that  the  concept  of  free  interaction  will 
yield  results  in  sufficiently  good  agreement  with  experiment  up 

to  the  reattachment  point  to  warrant  its  continued  development. 

In  any  event,  it  can  serve  as  the  zero-order  approximation  for  a 
more  precise  theory  which  takes  fuller  account  of  upstream  influ¬ 
ence  and  the  basic  elliptic  nature  of  the  problem. 

(13)  It  was  found  necessary  to  match  the  second  derivatives 
of  the  inner  and  outer  velocity  profiles  at  their  point  of  junc¬ 
ture  to  obtain  stable  calculated  solutions. 
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10.  RECOMMENDATIONS  FOR  FUTURE  WORK 

It  is  believed  that  the  present  work  is  capable  of  broad 
extension  and  exploitation.  A  number  of  specific  suggestions  for 
future  work  follows: 

(1)  The  calculative  program  should  be  extended  to  axially 
symmetric  flow  under  the  assumption  that  the  boundary  layer  is 
thin  compared  to  the  radius.  This  the  authors  intend  to  do. 

(2)  The  analysis  is  capable  of  extension  to  nonadiabatic 
surfaces.  Such  an  extension  is  contemplated  as  continuing  work 
under  the  present  program. 

(3)  Higher  approximations  within  the  framework  of  the  present 
theory  should  be  attempted.  These  approximations  should  consider 
additional  moments  of  the  momentum  relationship  for  the  inner  flow 
as  well  as  continuity  of  velocity  derivatives  higher  than  the  sec- 
on1  at  the  boundary  between  the  inner  and  outer  flows. 

(4)  Extension  of  the  calculative  method  downstream  of  the 
reattachment  point  is  desirable.  For  this  purpose  it  appears 
desirable  to  rotate  the  coordinate  system  perpendicular  to  the 
ramp. 

(5)  An  investigation  should  be  made  to  extend  the  analysis 
to  higher  ramp  angles. 

(6)  The  computational  program  should  be  used  as  an  analog 
device  to  study  separated  flows  systematically.  Effects  of  Mach 
number  and  Reynolds  number  should  be  studied  as  well  as  reattach¬ 
ment  flow  and  incipient  separation. 

(7)  A  systematic  comparison  between  the  results  of  the 
present  calculative  method  should  be  made  with  all  available 
exper imental  data  to  assess  its  limitations  and  the  necessity  of 
higher-order  approximations . 

(8)  The  present  method  should  be  capable  of  simple  exten¬ 
sion  to  blowing  or  suction  in  laminar  boundary  layers  including 
the  separated  region. 

(9)  It  would  be  of  interest  to  apply  the  method  to  turbulent 
boundary  layers  through  the  use  of  the  eddy  viscosity  approach. 
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(10)  The  present  results  for  a  Prandtl  number  of  unity  should 
be  generalized  to  other  Prandtl  numbers  and  to  real  gases. 

(11)  Since  the  Blasius  profiles  used  as  initial  conditions 
may  be  unrealistic  at  high  Mach  numbers  or  low  Reynolds  numbers 
because  of  induced  pressure  effects,  an  attempt  should  be  made  to 
start  the  calculative  program  at  the  leading  edge. 

(12)  Well-designed  careful  experiments  are  badly  needed  to 
obtain  detailed  flow  information  in  the  separated  region  at  both 
separation  and  reattachment  for  evaluating  theory. 

(13)  As  a  result  of  the  numerical  calculations  of  flow  in 
the  neighborhood  of  the  reattachment  point,  certain  speculations 
concerning  a  reattachment  criterion  have  been  made.  Further  work 
is  required  to  develop  a  sound  reattachment  criterion  both  from 
the  theoretical  and  experimental  points  of  view.  Theoretical 
work  is  needed  for  wedges  of  finite  length  and  infinite  length. 
Experimental  studies  of  the  flow  in  the  region  of  reattachment 
are  definitely  needed  to  guide  the  theory.  Accurate  velocity 
profiles  are  needed,  as  well  as  measurement  of  reattachment-point 
locations . 
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TABLE  II 

COMPARISON  OF  BLASIUS  PROFILE  WITH 
DORODNITSYN  APPROXIMATION 


_a_ 

VI 

Blasius 

Dorodnitsyn 

u 

VT  du/dT) 

u 

O.C 

0.0000 

0.3321 

0.0000 

0.3322 

.2 

.0664 

.3320 

.0665 

.3330 

.4 

.1328 

.3315 

.1331 

.3325 

.  6 

.1989 

.3301 

.  1  >94 

.3305 

.8 

.2647 

.3274 

.2652 

.3271 

1.0 

.3298 

.3230 

,3302 

.3220 

1.2 

.3938 

.3166 

.3939 

.3151 

1.4 

.4563 

.3079 

.4561 

.3064 

1.6 

.5168 

.2967 

.5163 

.2956 

1.8 

.5748 

.2829 

.5742 

.2826 

2.0 

.6298 

.2668 

.6292 

.2675 

2.2 

.6813 

.2484 

.6810 

.2501 

2.4 

.7290 

.2281 

.7291 

.2307 

2.6 

.7725 

.2065 

.7731 

.2094 

2.8 

.8115 

.1840 

.8128 

.1867 

?,0 

.8461 

.1614 

.8478 

.1632 

3.2 

.8761 

.1391 

.8781 

.1398 

3.4 

.9018 

.1179 

.9037 

.1172 

3.6 

.9233 

.0981 

.9250 

.0961 

3.8 

.9411 

.0801 

.942  3 

.0773 

4.0 

.9555 

.0642 

.9561 

.0610 

4.2 

.9670 

.0505 

.9669 

.0474 

4.4 

.9759 

.0390 

.9752 

.0363 

4.6 

.9827 

.0295 

.9816 

.0275 

4.8 

.9878 

.0219 

.9863 

.0206 

5.0 

.9916 

.0159 

.9899 

.0154 

5.2 

.9943 

.0113 

.9926 

.0114 

5.4 

.9962 

.0079 

.9946 

.0084 

5.6 

.9975 

.0054 

.9960 

.0062 

5.8 

.9984 

.0037 

.9971 

.0046 

6.0 

.9990 

.0024 

.9979 

.0033 

6.2 

.9994 

.0016 

.9984 

.0024 

6.4 

.9996 

.0010 

.9989 

.0018 

6.6 

.9998 

.0006 

.9992 

.0013 

6.8 

.9999 

.0004 

.9994 

.0010 

7.0 

.9999 

.0002 

.9996 

.0007 

7.2 

1.0000 

.0001 

.9997 

.0005 

7.4 

1.0000 

.0001 

.9998 

.0004 

7.6 

1.0000 

.0000 

.9998 

.0003 

7.8 

1.0000 

.0000 

.9999 

.0002 

8.0 

1.0000 

.0000 

.9999 

.0001 

8.2 

1.0000 

.0000 

.9999 

.0001 

8.4 

1.0000 

.0000 

1.0000 

.0001 

8.6 

1.0000 

.0000 

1.0000 

.0001 

8.8 

1.0000 

.0000 

1.0000 

.0000 
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TABLE  III 

VALUES  OF  C0,  c1,  C2,  AND  c  FOR  CALCULATIVE 

3 

EXAMPLE  CORRESPONDING  TO  LEIGH  CASE 


X/i 

Co 

ci 

C2 

°3 

0.0100 

0.0100 

0.3157 

-0 

.1923 

-0 

.0313 

1.1000 

.0403 

.0404 

.6174 

- 

.3806 

- 

.0565 

.9991 

.1005 

.1011 

.9350 

- 

.5829 

- 

.0775 

.8550 

.1485 

.1499 

1.0980 

- 

.6881 

- 

.0854 

.7495 

.2022 

.2049 

1.2325 

- 

.7744 

- 

.0906 

.6425 

.2509 

.2549 

1.3255 

- 

.8328 

- 

.0939 

.5559 

.2969 

.3027 

1.3959 

- 

.8753 

- 

.0967 

.4821 

.3507 

.3589 

1.4616 

- 

.9125 

- 

.0998 

.4055 

.4019 

.4125 

1.5117 

- 

.9381 

- 

.1035 

.3409 

.4480 

.4613 

1.5488 

- 

.9543 

- 

.1073 

.2890 

.5094 

.5268 

1.5889 

- 

.9677 

- 

.1136 

.2281 

.5555 

.5763 

1.6134 

- 

.9723 

- 

.1193 

.1878 

.6016 

.6261 

1.6341 

- 

.9729 

- 

.1261 

.1517 

.6528 

.6818 

1.6535 

- 

.9693 

- 

.1351 

.1159 

.  7040 

.7380 

1.6699 

- 

.9618 

- 

.1456 

.0845 

.7501 

.7890 

1.6828 

- 

.9520 

- 

.1567 

.0596 

.8013 

.8460 

1.6960 

- 

.9386 

- 

.1708 

.0356 

.8525 

.9035 

1.7095 

- 

.9242 

- 

.1863 

.0155 

.8986 

.9556 

1.7265 

- 

.9182 

- 

.1978 

.0014 

.9011 

.9585 

1.7279 

- 

.9189 

- 

.1978 

.0008 

.9024 

.9600 

1.7288 

- 

.9194 

- 

.1977 

.0005 

.9037 

.9614 

1.7297 

.9202 

.1976 

.0002 
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TABLE  IV 

COMPARISON  OF  LEIGH  PROFILE 
WITH  PRESENT  SOLUTION 

(a)  X^/i  -  X/i  -  0.000142 


U/U 

o 

— 

Leigh 

Dorodnitsyn 

0.0 

0.0000 

0.0000 

.2 

.0026 

.0027 

.4 

.0096 

.0099 

.6 

.0210 

.0216 

.8 

.0366 

.0376 

1.0 

.0564 

.0578 

1.2 

.0802 

.0819 

1.4 

.1078 

.1098 

1.6 

.1390 

.1412 

1.8 

.1735 

.1757 

2.0 

.2110 

.2130 

2.2 

.2510 

.2527 

2.4 

.2931 

.2944 

2.6 

.3368 

.3376 

2.8 

.3815 

.3819 

3.0 

.4266 

.4267 

3.2 

.4716 

.4714 

3.4 

.5158 

.5156 

3.6 

.5586 

.5586 

3.8 

.5995 

.5999 

4.0 

.6381 

.6389 

4.4 

.7065 

.7082  ! 

4.8 

.7619 

.7640 

5.2 

.8039 

.8054 

5.6 

.8338 

.8340 

6.0 

.8536 

.8525 

6.4 

.8658 

.8638 

6.8 

.8729 

. .8706 

7.2 

.8767 

.8745 

7.6 

.8787 

.8767 

8.0 

.8802 

.8780 

Dorodnitsvn 

X  /i  -  0.962814 
8 

X/i  -  0.962672 

c  -  1.730535 
o 

cx  -  -0.92109 
c 2  -  -0.197233 
c3  -  2.32X10"5 
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TABLE  IV.-  CONCLUDED 

(b)  X  /i  -  X/i  -  0.002542 
s 


U/U 

o 

Leigh 

Dorodnitsyn 

0.0 

0.0000 

0.0000 

.2 

.0041 

.0042 

.4 

.0126 

.0129 

.6 

.0255 

.0260 

.8 

.0426 

.0435 

1.0 

.0638 

.0650 

1.2 

.0890 

.0905 

1.4 

.1180 

.1197 

1.6 

.1504 

.1522 

1.8 

.1861 

.1877 

2.0 

.2246 

.2260 

2.2 

.2655 

.2665 

2.4 

.3083 

.3089 

2.6 

.3524 

.3526 

2.8 

.3974 

.3972 

3.0 

.4427 

.4421 

3.2 

.4875 

.4868 

3.4 

.5314 

.  5307 

3.6 

.5737 

.5732 

3.8 

.6138 

.6138 

4.0 

.6515 

.6520 

4.4 

.7178 

.7193 

4.8 

.7708 

.7726 

5.2 

.8106 

.8116 

5.6 

.8384 

.8382 

6.0 

.8566 

.8552 

6.4 

.8678 

.8656 

6.8 

.8741 

.8717 

7.2 

.8775 

.8753 

7.6 

.8792 

.8773 

8.0 

.8805 

.8785 

Dorodnitsvn 

X  /l  -  0.962814 

s 

X/i  *  0.960272 

c  -  1.728945 
o 

cx  -  -0.919618 
c2  -  -0.197710 
c 3  -  4.588X10"4 
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TABLE  V (a) 


COMPARISON  OF  DORODNITSYN,  QUARTIC,  AND 
SIMILARITY  VELOCITY  PROFILES 
AT  SEPARATION 


Y/6  * 

Dorodnitsyn 

u/u1 

Quartic 

u/ui 

Similarity 

u/ui 

0 

0 

0 

0 

0.0848 

0.0046 

0.0066 

0.0040 

.1696 

.0183 

.0252 

.0158 

.3391 

.0719 

.0914 

.0636 

.5935 

.2098 

.2406 

.1927 

.8478 

.3961 

.4177 

.3802 

1.0174 

.5315 

.5368 

.5230 

1.2717 

.7226 

.7004 

.7278 

1.4413 

.8252 

.7927 

.8364 

1.7804 

.9470 

.9252 

.9598 

1.9500 

.9736 

.9644 

.9839 

2.6282 

.9987 

1.0000 

.9999 
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TABLE  V (b) 

COMPARISON  OF  DORODNITSYN  AND  QUARTIC  VELOCITY 
PROFILE  DERIVATIVES  THROUGH  THE 
BOUNDARY  LAYER  AT  SEPARATION 


Y/6* 

Dorodnitsyn 
d (U/Ux ) 
a (Y/6*) 

Quartic 
d  (U/Ux ) 
d (Y/6*) 

0 

0 

0 

0.0848 

0.1081 

0.1519 

.1696 

.2144 

.2829 

.3391 

.4141 

.4865 

.5935 

.6554 

.6627 

.8478 

.7882 

.7110 

1.0174 

.7977 

.6870 

1.2717 

.6786 

.5894 

1.4413 

.5256 

.4962 

1.7804 

.2095 

.2832 

1.9500 
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FIGURE  I.  -  SEPARATED  LAMINAR  SUPERSONIC  FLOW  OVER  A  FLAT  PLATE 
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PRE-SEPARATION  REGION. 
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FIGURE  5.  -  PRE-SEPARATION  PRESSURE  DISTRIBUTIONS  FOR  INITIAL 
PRESSURE  DISTURBANCES  OF  VARIOUS  MAGNITUDES. 
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FIGURE  6.  -  VARIATION  OF  Cn,  C,,  C,  ,  ANO  C,  WITH  AXIAL  DISTANCE. 
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FIGURE  7  -  EXAMPLE  VELOCITY  PROFILES  GIVEN  BY  EQUATION  (96). 
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FIGURE  8.  -  BEHAVIOR  OF  Co,  C|,  C2,  C3,  GOING  THROUGH  SEPARATION  POINT 
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FIGURE  9.  -  SKETCHES. 
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(o)  VELOCITY  PROFILES. 

FIGURE  II.-  VELOCITY  PARAMETERS 
FAR  DOWNSTREAM  OF  SEPARATION 
ON  BASIS  OF  VARIOUS  ASSUMPTIONS. 
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FIGURE  13-  NONDIME  NSIONAL  VELOCITY  PROFILES  DEVELOPED  ASYMPTOTICALLY  BY  SOLUTION 

FOR  CASE  C. 
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FIGURE  14.-  CALCULATED  FLOW  PATTERN  AND  PRESSURE  DISTRIBUTION  IN  REATTACHMENT  REGION. 


CASE  C 


FIGURE  15  -VELOCITY  DISTRIBUTION  THROUGH  THE  BOUNDARY  LAYER  IN  THE  REGION  OF 

THE  COMPRESSION  CORNER  FOR  CASE  C. 
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EFFECT  OF  FAIRING  RADIUS  ON  CALCULATED  PRESSURE  DISTRIBUTIONS. 


FIGURE  18.- EFFECT  OF  INNER  VELOCITY  PROFILE  APPROXIMATION  ON  CALCULATED  PRESSURE 

DISTRIBUTION. 
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FIGURE  19.- CONTINUED. 


2.0 


fEEN  MEASURED  AND  CALCULATED  PRESSURE  DISTRIBUTIONS 
FOR  FLAT  PLATE  WITH  10°  WEDGE. 
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FIGURE  20  -  CONCLUDED 


Ro/  2  =  8.5  X  lOVFT. 


IN  BEGINNING  OF  INTERACTION. 
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FIGURE  22  -  VARIATION  WITH  DOWNSTREAM  DISTANCE  OF  OSWATITSCH 
REATTACHMENT  ANGLE  AND  DIVIDIN6  STREAMLINE  SLOPE  (RELA¬ 
TIVE  TO  WEDGE). 
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FIGURE  23-  VARIATION  OF  OSWATITSCH  REATTACHMENT  ANGLE 
FOR  DIFFERENT  POSITIONS  OF  BEGINNING  OF  INTERACTION. 
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FIGURE  24.-  VARIATION  OF  QUANTITIES  ALONG  DIVIDING  STREAMLINE  APPROACHING 

REATTACHMENT. 
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(b)  MACH  NUMBER  ON  DIVIDING  STREAMLINE. 


AFFDL-TR-65-107 


APPENDIX  I 

CALCULATION  OF  g  FUNCTIONS 


The  family  of  functions  of  a  single  variable  occurring  in 
the  analysis  are  given  by  Equation  (43)  as  follows: 


g  (c  )  =  u  -qu-  ,  c  >  0  (1-1) 

n  3  J  ~\l  u  +  c 

O  V  3 


Simple  integration  yields, 

90(C3)  =  2(VX  +  S  ‘  VS)  <I_2) 

The  other  g functions  can  then  be  simply  generated  by  the 
following  recursion  relationship 


9 


n 


2 

2n  +  1 


(V1  +  S 


nc 


3 


g 


(1-3) 


For  large  values  of  c3  the  recursion  formula  requires  many 
more  decimals  of  calculation  than  significant  figures  in  the  cal¬ 
culated  result  so  that  a  series  solution  rapidly  convergent  for 
large  c3  is  desirable. 

The  following  results  are  used  for  large  c3 


g  (c  ) 
n  3 


i 


o 


n  , 
u  du 


du 


(1-4) 


g  (c  ) 

^n  3 


1 

n  +  1 


+ 


jk 

n  +  k  +  1 


k=  l 


(1-5) 


=  (-pk(2k  -  i): 
22k_1kl (k  -  1)  : 


(1-6) 
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APPENDIX  II 

DERIVATION  OF  FREE -INTERACT ION  RELATIONSHIP 


The  basic  assumption  for  calculating  the  pressure  gradients 
of  the  external  flow  field  is  that  such  gradients  are  the  same  as 
those  for  an  inviscid  flow  with  the  solid  boundaries  augmented  by 
an  amount  5*  due  to  the  boundary  layer.  The  particular  invisid 
flow  relationship  used  for  the  calculation  is  not  germane  to  the 
basic  assumption,  and  in  this  case  was  taken  to  be  the  Prandtl- 
Meyer  shock-expansion  relationship 


P 


1 


yM. 


COS  0 


d  0 


cos20 


(II-l) 


It  is  noted  that  in  free  interaction  the  pressure  distribution  is 
thus  dependent  on  the  variation  of  6*  along  the  body,  and  the 
variation  of  6*  in  turn  depends  on  the  pressure  distribution. 
This  interrelationship  will  be  termed  the  free-interaction  equa¬ 
tion.  Specifically,  it  will  consist  of  two  equations,  Equation 
(II-l)  and  the  following  equation  relating  0  to  the  boundary- 
layer  quantities. 


db  *  db  .  * 

o  ,  i  ,  dw  _ 

dx  dx  dx 


tan 


0 


where 

bQ*  displacer  ant  thickness  of  outer  flow 
displacement  thickness  of  inner  flow 
dw 

slope  of  streamwise  boundary 


(II-2 ) 


Pre-Separation  Region 

The  calculation  of  b  *  and  its  derivative  is  a  fairly 

o  1 

laborious  algebraic  operation  which  will  now  be  carried  out.  The 

definition  of  6  *  is 

o 
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00 


We  will  transform  the  integral  to  t. 


00 


Now  since 


1 


u  =  U 


a 

o 


we  have 


(1  +  mi) 


(1  +  )  -  m^u2 


so  that 


5  * 

o 


P0ao 

Plal 


[<! 


m1 )  (1  - 


u)  + 


mxu  (1 


-  u)' 


In  the  outer  flow  the  velocity  profile  is  given  by 


(1  -  u)  y  u  +  c3 
C0  +  C!U  +  C2U? 


(II-3) 


(II-4) 


(II-5) 


dr,  (II-6) 
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with  the  result  that 


rU  t v 


5  *  =  — °-Q  — 77^ — —  (1  +  m  )  (c  g  +  c_  g  +  c  g  )  +  m.  (c  g_  +  c  g_  +  c„g  ) 

O  O  a  U  I  1  O  O  11  2^2  1  O  1  12  233 

'll  i  -1 


(II-7) 


We  now  form  the  quantity 


Rol/2  d5o*  ,  /  1  di  dX  X  f(1  + 

1  +  mQ  dx  o  l  1  +  mQ  J  dX  dx  d£  UJL 


'(1  +  mi^Cogo  +  Cigi  +  C2g2^  +  mi  (Cogi  +  Cig2  +  C2g3^ 


(II-8) 


Carrying  out  a  tedious  differentiation  and  noting  that 


dmi  .  ,  ui 

3T  *  "V  =  2l"i  5; 


(IX-9) 


yields 

Rol/2  d6Q*  /  n>i  \  /  mi  \  .  f  mi  \ 

1  +  mQ  dx  ~  (9o  +  gi  1  +  mx  J  co  +  (^1  +  g2  l  +  m^y0!  +  (^2  +  g3  1  +  J  ^2 

mi 

+  cg+cq+cg)  +  -r- ; -  (eg  +cg  +cg)  c 

0^0  1^1  2^2'  1  +  m  oyi  1^2  2^3  3 

1 

,  mi  f.  ,  ,  .  At  -  1  1  +  mi  \ 

+  r~T^  <cogo  +  cig!  +  c^2)  r r^T  '  ) 


+  (c  n  +  c  g  +  c  g  >  1  11 

o  1  12  2^3  y  -  1  1  +  J 

„  1/2  (tan  *  -  fr) 


1+m 


(11-10) 


Equation  (11-10)  taken  together  with  Equation  (II--1)  constitutes 
the  free-interaction  relationship. 
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Post-Separation  Region 


After  the  separation  point  the  inner  flow  causes  the  turning 
of  the  inviscid  flow  by  an  amount  depending  on  6^*.  By  analyzing 
Equation  (II-6)  we  can  viite  an  equation  for  5^*  as  follows: 


p  a  I  r 

6-*  =  io_o  _y — a — o  r( 

1  Plal  U1  J  L 


1  +  m^Hl  -  u)  +  miu(l  -  U)]  dr)  (11-11) 


In  this  case  the  velocity  profile  is  given  by 


/a  +  2a  \  /a  +  a  \ 

—  /  s  W  \  2  ,  /  S  W\3 


(II-12) 


The  resulting  equation  for  5^*  is 

v  ^  _l_  /i  ^1yi/2(T-i) 

1  ~  Ul  R  ^  V1  +  mo  J 

o 


(!  +  ml)  rj  - 


i r  (aw  -  as’ 


/  a  2  a  a  a  2 

n  3  _w_  _  JtLJL  +  _s_ 

s  V  105  70  105 


where 


£u  V  d 

a  =  =  - 

s  oti  c 

^s  ° 


(11-13) 


(11-14) 


Carrying  out  a  differentiation  as  before  yields 
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Introduction  of  db^*/dx  in  Equation  (II-2)  along  with  d&Q*/dx 
yields  the  f ree-interaction  equation  when  taken  together  with  Equa¬ 
tion  (II-l)  .  Note  that  the  wall  may  have  any  arbitrary  shape  given 
by  w(x).  It  case  a  quadratic  rather  than  a  cubic  is  used  for  the 
inner  velocity  profile,  we  set 


a 


w 


a 

s 


Vs. 

c 

o 


(11-16) 
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